UNCLASSIFIED 


AD NUMBER 
NEW LIMITATION CHANGE 
TO 


Approved for public release, distribution 
unlimited 


FROM 


Distribution authorized to U.S. Gov’t. 
agencies and their contractors; 
Administrative/Operational Use; Feb 1961. 
Other requests shall be referred to Office 
of Naval Research, 800 North Quincy 
Street, Arlington, VA 22217-5660. 


AUTHORITY 


ONR ltr, 15 Jun 1977 


THIS PAGE IS UNCLASSIFIED 


; : OAS: IFT TED _ 


he REN 


gy 254 412 © 


by the 


» 7 HERVICES TECHNICAL INFORMATION AGE 
ARLINGTON HALL STATION 
ARLINGTON 12, VIRGINIA 


UNCLASSIFIED 


e 


NOTICE vermment or other dravings. specl- 
Cicaiion co csther data are used for ::." ;urpose ' 
other <hca, connection with a «.finite.y related 
governns ¢  -ocurement operation. the U. &. 
Governme t 1 ereby incurs no rerficsib!* {t;, nor ary 
obligat’ on. wztsoever; and the 242. thal ua? Gover- 
mert. may .av formulated, furnishea, or in any way 
suppilec che said drawings, specifications, oi other 
data {s ot t) be regarded by implication ->- sther- 
Wie. 7 no© ° Magner Licensing the oolde: .7 any 
wa 8 * OF COrporati ua, wr comnvering.@ , rights 
or perm: s.on to manufacwure, use or sell ezy . 
“at te.) invention that may in anv way be releted 
the eto. 


aot , 
a ence rc a a a 
- = = To a = = eee r ene = : a emer Al’ ee . - J. - IF. 
ate = 7 — 7 ike a Paha See, cule eter 
cast ea : ’ e “re Dolan 
ee u es Pan Y =; a 4 igs ~~ 
: | : ee ee “8 
k toy an = ese , : 
NS EN ea nar PES Exe ES ete Scie és ee at ) : 1 re 
. _ x ; 


es 


lulerim Report 


T-A2049-12 
4) AN ANALYTICAL AND EXPERIMENTAL 
ie STUDY OF THE STABILITY OF EXTERNAULY- 
A. j PRESSURIZED, GAS-LUBRICATED 
Z +. i THRUST BEARINGS 
a 
re | 
<Z. J by 
x. L. Licht and H. G. Elrod, Jr. 
rail eee 1 
Css February 1961 a 
rare bee een 
CD ene 
Ty ie Prepared under 
=L ec 
— ad Contract Nonr-2342(00) 
<o ee Task NK 061 - 113 rae a oe Oe, 
wed BT toa tf 
Jointly Supported by . 
DEPARTMENT OF DFFENSE oe 
ATOMIC ENERGY COMMISSION 
MARITIME ADMINISTRATION ae 
NAYEONAL ARRONALUTUCS AND SPACE ADMINISTRATION | P i 
Administered by 
OFFICE OF NAVAL RESEARCH 
Department of the Navy Copy No. 


oF % : : , i me. Wo . 


~~ « ae r 


, 26 / Bee 
THE FRANKLir iNSTITUT 
LABORATORIES FOR RESEARC'. AND DEVELOPME 
PHILADELPHIA ; PENNSYELVAP 


THE FRANKLIN iNSTITUTE .- Laborarsiies for Research and Develot:ment 


Interim Seport I--A2049- 12 


AN ANALYTICAL AND EXPERIMENTAL 
STUBY OF TRE STABILITY GF EXTERNALLY- 
PRESSURIZED, GAS-LUBRICATED 
THRUST BEARINGS 


by 


I. Licht and H. G. Elrod, Jr. 


Prepared under 


Contract Nonr 2342100) 
Task NR 061-113 


Jointly Supported hy 


DEPARTMENT OF DEFENSE 
ATOMIC LNERGY COMMISSION 
MARITIME ADMENISTRATION 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Administered by 


OFFICE OF NAVAL RESEAKCH 
Department of the Navy 


Reproduction in Whole or in Part ts Permitted 
for any Purpose of The U S_ Government 


THE FRANKLIN INSTITUTE © Laboratories for Research and Development 


T-A2049- 12 


FOREWORD 


This ceport deals with both the thenretical and experimental phasca of an iavesiigation of 
the stability of externally-pressurized, gas-lubricated bearings. 


The analysis, presenied in Interim Report No. 1-42049-4, has been extended and comple- 
meated by additional material in Appendices 2 through 7 and i» included in this rport for the 
sake of completeness. Detailed descriptions of apparatus, instrumenis and experinz:+ ure give 
and are supplemented by €5 figures, containing graphs drawings and photographic records. ' 


The theoretical predictions of the present and of earlier analys“s are compared with ex- 
perimental results, areas of further resecich are suggested and recommendations are made for 
utilizing the results of the present investigation for desiyn purposes, 
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1, ABSTRACT 


The subject of this investigation is the stability of externally-pressurized, gas-lubricated 
bearings. An analysis is made which cats the flow within the gas film on a distributed, rather 
than on a lumped parameter basis, Rectangular and circular bearing configurations are considered 
and the effect of several parameters on stability is investigated. Results uf the present analysis 
are compared .ind contrasted with those obtained from earlier, simplified theories. The experi- 
mental program pursue. here is guided by the analvsie. An apparatus ia designed, constructed 
and instrumented to provide for the variation of goveming parameters and for the accurate deter- 
mination of their critical values. Comparison of theory and experiment is made and the validity 
of theoretical assumptions, relative to the experimental model, is discussed. 
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He INTRODUCTION 


Instubiitties af Gast ubricated Rearings and Related Problems — Review of 


* Pertinent P iterature 


Reewnt abv osees iy space, computer and nuclear technology have precipitated a renewed 
interest in pas hibeication, Gas bearings are being successfully incerporated in systems pre- 
cicding the use ef cons. mtivnal lubricants and rolling elements Vo quote contrasting examples 
of application ~ ges bearnas have been used for small turbine rotors, driven and supported by 


False at very low temperatures in gas liquefaction machines. ‘They have 
hod of gveoscope support in inertial guidance devices, a8 sensing elements 


he toesxhaust Leta og ad 


been employed ase 
ia computers and as lowefriction guides in accelerometers, dvaamometers, vibration tables and 
torque sieasuring insivuments, Ip nuclear reocturs, gas bearings have provided a solution for 

sealed and isolated +: stems, in whieh radioactive gases furnish the only available lubrication 


re 
THe a yay 
Gas bearings mav be broadiv classitied into two types: 


a. Selfeceting bearings. in which the load capacity is induccd through the retative motion 


uf suefaces, 


bo Uscernalls-pressurized bearings, which rely for their load capacity or an external 


pressure Source, 


One 


aive difficulty encountered in applications of gas bearings is their susceptibility to 
selfeexcited vibrations. In unstable, self-acting bearings, the center of the mtating journal may 
cuher describe a closed path within the clearance ciecle, or follow a spiral path until surface 
contact and destruction take place. This type of instability, generally ceferred to in literature as 
“whirl, occers also with incompressible lubricants, but appears to be particularly troublesome 
with gases. It han been the subject uf numerous theoretical and experimental investigations 

125, 37, 38. 90, 10, 11, 42, 4317, but its mechanism is not yet fully understood. 


The subject of the present investigation is un instabitity, associated with externally- 
pressurized bearings and sometimes referred to as ‘‘air’’ or “‘pneumatic hammer’. It is fre. 
quently encountered in valves, regulators and other elements of nreumatic and hydraulic contro! i 
systems (8,9, 12, 13, 14 15, 16], and is closelv allied to such seemingly unrelated phenomena 
as self-excited vibrations of fuel injecters (17, 57] and pressure reducers [52], axial auto-oscil- 

. lations of steam-cushioned turbine rotors [57! and vinlent pressure flucuations in fan-ventilated 
chambers {18 57]. The ‘pneumatic hammer’’ in externally-pressucized bearings is not associated 
with the steady, relative motion of surfaces, nor greatly influenced by it at moderate rotational 
speeds. If, however, the amplitude of vibration becomes excessively large in relation to the 
clearance, surfaces mtv contact and scoring, or complete destruction ensue. Vorenver, even 
relatively ninute auto-oseillations cannot be tolerated in applications involving sensitive instru- 
nents, or inertial guidance devices, On the other hand, unless the selection of geometrical and 


a a 


Tt Numbers in bractets refer to bibliography. i 
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ea : bo kek ; : ae 
fluid oarurctors rs governed $ at beast > qualitative understanding of the physical mechanism of 


wnstabilitn dt is rrcheble that the acing with be unstable in ail, or part of the contemplated 


region of operation. 


The phetomenan was observed in externally-pressurized thrust and journal bearings by 

Po Caaetet 1p, Pigott and Vucks (238 Fischer, Chetshim and Fullec 
aed Sterns PSL Graneck and Kere £304, Ailea, Stokes and Wiitley [321 and many 
Prygt establish a ; 


sacctaal guidwee. A unmber of investizators have speculated on likely causes of instability and 


ufhees 


urical criteria for specific bearing configurations used in 


sstec renndial neasuees. based on ohsecrvation, of on intuition 1%, 30, 32 S02 The 


first. rational exnlanations of the phenamenon in gas lnbrication were given independently by 


have s 


th Palier and Stend:cht $51. whe considered the instability of a simple, air-lubricated thrust 
beaing, and by Richardson [2. 51, who concentrated on the dynamics of externaliy-pressurized, 


nw the initial eccentricity te he small and motion tu be 
More recent!v. Licht [T! extended the analysis to multiple- 
journal bearings, assuming plane mocon of the journal, initially in an eccentric equilibrium 
posttien, Roudebush fé] considered the stabslity of a circular thrust bearing and obtained digital 


aie-lubric ated journal bearings. assn 


simited ts a single degree of freedom. 


computer solutions for the nonlinear, lumped, fluid-dvnamic equatioas. [fis conclusions were in 
agreement with “hose ‘smsed on linearized analyses, altL ough compater solutions were tecessarily 


fimited to a oumver of specific canes. 


the analyses referred to in the preceding paragraphs were of the lumped-parameter type 
aud celied on arbiteary, a priori assumptions with recard to variation of nreseure profile and flow 
in the bearing clearance. The objective in making these assumptions was to permit a simple 
formulation of the problem and te avaic the use of partial differential equations in the mathe- 
matical description of the fluid-dvnamic phenomena occurring in the bearing clearance. 


2. feepose and Outline of the Present Investigation 


Because of the large number of parameters which influence the stability of extemnally- 
pressurized bearings, and because of the complexity introduced through nonlinearity, the present 
investiga 


nis restricted to a perturbation analysis. The limitation that perturbation analyses 
give necessary conditions of instability, not suificient conditions of stability, is recognized, 
Lut linearization hi» frequently provided a good first anpreximation in the solution of fundamente! 
problems of fluid dynamics. More specifically, analog and digital computer studies of pneumatic 
controls indicate that results obtained from linearized equations often compare very favorably 
with exact solutions (10, 52]. Some encouragement may also be derived from the fact that analog 
solutions of non-linear equations for a Leariag using an incompressible Ibricant, but geometri- 
cally identical with the mode! pmpesed in this investigation. were in good agreement with results 
obtained from their linearized «sui ents FLL]. 


Results vf previous experiments, reparted in reference [4], agreed only qualitetively with 
theoretical predictions of simplified analyses [4, 5]. It was felt that, although the sinplified 
analyses explain. a the basic aspects of the physical mechanism of ‘‘pneumaric hummer’ is gas 


hearings the dearna nf aimnlifinatinn in fgrmnlation the orobt 
bits : ; ; qrme 


SSE ee praaacs 


ee eee a ee a ee 
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[2, 3] recognized that dynamic changes in pressure-profile depended on time rates of boundary 
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vation. tn fact, # priori assumptions with regard to pressure-profiie variation and mass contin- 
wity appear to have little justificatise, particularly for corfigurations in which the surface cor- 
responding to the slot coustitutes a large fraction of the tota! bearing area. Equally questioa- 
ubte are the results of simplified ana!vses, according to which stability is not influenced by the 
hearing mass 13, 51.f 
' 
In the theprgatical part of the present icvestigation, while retaining the method of small 
petturbations vad'a cuuber of simpiifving assumpiions, the nousteady flow ix the bearing clear- 
“orce is treated on a conlingous, rather than on a lumped basis. The fluid-dvnamic phenomenon . 
sin the slot is descrihed by a partial differential equation. The latter is required to satisfy the 
boundary condition at the outer bearing perimeter and to be compatible with the continuity 
equation for a centrally facated, variable-volume recess, and also with the equation of motion of 
the beaeing asians. Stability analyses are made for axial motion of parallel-surface thrust bearings. 
Rectangular and circular conficurations’t ere considered and stability boundaries are obtuined to 
show the effect of several, important fearing parameters and to shuw luw ue results compare 


with those obtained froin lumped-paraneter theories. 


Ja the expermenta’ part of the present inveetiaation, the bearing is circular and identical 
in geometry to the theoretical model, The enrface bounded by the perimeter of a centrally- 
located, adjustable-wolume recess, constitutes a saul! fraction of ihe total bearing area, in order. 
to hinder the oscene of fluid and magnify the trapping and squeezing effect in the slot, Provision 
is also mide to vary the bearing mass without altering the lead, in order to observe the influence 
of mass on stability. 


een ce ee Re ctr 


wee a er a a 


T See also Appendix 2. 


tt Schematics, diawings and photographs of actual thruat plates are shuwg in Figures I, 7, 8, 13, 19 ; ‘ : 
and 20. i 
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It. aNatysis? 


Horivation of (ianracteriatic Pa stiohs 


fii chat folicws. equations are simlianeansty devolone i for the rectangular and circular 


izue ations stewn in} igures Leand ib, Fquations applicable to both geometries are denote ! 
; vets. Those denoted by the letter a"? refer to the rectangular geometry, while 


teeoe donted by the letter ‘b'' refer te the crecular geometry. 


Paitations C1) throagh (1 ace those developed by G.L Taylor and PF. G. Saffman [1]. who 
Investisated comoressthility effects at iow Revnolds number, One case they considered consisted 
of tive parallel discs. separated Sy a parrow aie gar. when rotation was accompanied by small 
Mitude. forced vibsation in the axia: direction. 8v contrast, the present investigation cencerns 
itself with the onset of self-evcited oscillations. tt differs alsn becanse, in the abscuce of rota- 


tion tiuid flow is induced by means of an external pressure source connected to an internal, vari- 


ableevelame reservoir. 


it is assumed that forces due to inertia are smali compared with those due to viscous 
Stresses and that the pressure gradient i. the teansvecse direction of the stot 13 negligible. The 
weapditiele of oscillacions in the anticipated frequency band width can be considered arbitrarily 
smatlia cespartson with the magnitude of the gap widih. ‘Se mean fiow veiocity in the annulus 


can then be obtained from the reduced Vavier-Stokes equation av weitten as: 


H? oP 
PN ak) 
i Igy 3x (la) 
aie ee (1b) 
2p dr 


Since the dimension of the annu'us sepurating the metallic surfaces ia small in comparison 
with the linear dimension of the slot. it is reasonable to assume that conductivity prevents any 
appreciable variztion in temperature and ihat the relation 


‘pb 
is = constant - (2) 
P 


tepresents the thermodynamic process to a sufficient degree of approxinution. 


= 

tMtatecial contained in this ection was presented at the 23rd National Conference o: the Arplied Vechanics 
Division, ASME, June 1960, and was published under the thle "A Study of the Stability of F xternally 

ti. Gud, Jey Journal of Appiied Mechanics, Vot. 27, %o. 2, 


Dewecurizad (uc De acing = 


June 1960, po. 250-258. © 


coe tee ee ee 


} mes 


THE FRANKLIN INSTITUTE « Laboratories for Research and Development 


T-A2049-12 


The continuity equation then becomes: 


‘ ie eo ~ Ba) 
UH ou) = 0 | 


Oia ee ha) = (24) 
et So or 


Substhution of equations (1) and (2) in the continsity equatian (3} yields a dynamic lubri- 


cation equation: 


in which F(t) «is the excitation function. 


1 7:.« effect of external damping in negligibly small 


mee 


2 (SEV oe SS 
Ox : (5-) y3 at wen : 
| Lo (ep aP ae 2 (pH) (4b) 
tdr\ ar H* ot 
. 
The baxie assumptions and the derivation of these equations are given in Appendix 6. 
F quation (4) must satisfy the boundary condition 
PCR, t) = Pay (5h) 
PCN, t) = Pap (Sb) 
and be compatible with the equation of motiont: 
PH 
m—— = 2X,(P,-Pro} + 2p (P-P,)dx + F(t) (6a) 
dr’ X, ; 
tal 2 - 
m de = wR, (P.-P,,) +2n (P=P,)rde + F(t) (6b) 
t ; Ry 


re a are cee AR npn 


“se 
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In addition, the continuity equation of the variable-velume recess is as follows: 
aM, 5 : 
ew AW, (7) 
dt mss 


Ia expressing dhe ctantities on the right hand side of equation (7), it is avoumed that the 
fiow can be treated as quasi-steady. Theoretical and experimental investigations of pulsating is 
fw theouch orifices were conducted by Schalt--Grunaw | 461, who discussed the conditions 
under which instantaneoas flow cates could be evaluated from steady-fluw equations and gave 
ctiteria for their applicability. “To express ihe left hand side of equation (7), it is necessary to 
describe an extremely compley process, associated with sinullaneods charging and discharging ot 
a satiable-volume container. In the present analysis it wil! be assumed that fluid properties 
threughout the reces# volume are uniform at aay instaat-of time and that the process may be ade- 
quately deccribed by the pressure-densitv relation P pY = constant. The quasi-static charging 
and discharging of an adiabatic container, for example, may be considered as isentropic to the 
first degree of approxmnation. perf. i miving af entering gas being assumed. Experiments con- 
ducted by Shearer {9} and Shinver ave wagner [561 also indicate that. for relatively rapid rates 
of charwing and discharging. the anaumptisn Pp = constant is reasonable. The effect of wave 
Propagation, heat tranafer and trrevercibilities. associated with the processes of charging and 
discharging of fixed and variable volume containers. cannot always be neglected [| 51.53. 54,55 [. 
However. assunptioas identical with those made in implementing equation 'T) have bern frequently 
and successfuily emptoved in related problems. in the field of nneumatic controls and regulators 
(9.14.42). Consequently, the quantities appearing in equation (7) are expressed as follows: 


y (9\t 
| ites z (6 =), | | (9a) : 
W, =- a (rp aan , . (9b) 


¥ She flow through the supply aozzle can be approximated clonely by @ quadrant of ua ellipse whenever 
F/R, 19 considerably atove the critical-presnare rétiv. Tuis in Fliegner’s approximation, and is con- 


venieat for computations! purposes. 


-1 
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Me = ON 84ND, . (10a) 
M, = aR? (5 +H), (10b) 
while the pressure-density relation for the recess volume is taken as: 
aed . 
P P. 
a -_— (11) : 
Peni Pro 
In the next step che perturbation quantities 
P-P, = p 
ros Paha. P, (12 
H-~H, =h 


und iheie derivatives are introduced. Neglecting their squares and cross products ae quaucities of 
higher order, we obtain instead of equations (4), (6) and (7) a set of corresponding expressions. 


‘The rontinear, partial differential equation (4) is replaced by one which is linear, but 
contains variable coefficients. 


Thus in terms of perturbation quantities we have: 
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!rom ejuation (6° we obtaia. 


in whieh Fret) as the disturbance, 


rquation (7) reduces to 


in which the coefficients hy. 
; 


iM, \ 
7 
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[ ak tv 
mo = 2X5 1 aA pdx + f(t) (14a) 
| i. wXy 
i a : | 
om d al = nite’ p, +is f iprdr + f(r) (14b) 
| dtr Up 
dh oe aa . OP, op, 
| h, = ~hyh + {Ry + Ks dp, +k, aoc he ra (15a) 
a 3 
j : dh . a a P, . 
tl C ao? Cyh = (C+ C)p, + Cx oo ~C, re (15b) 
and (. represent the quantities listed in the following: 
i rep q 8 
*% i 
“Min ee zs) soak 
ee ap 1, ve 
| xaXyp 
- a (ii 4 ') ky sims bs mee 2p 4, 
yFT, VH, 5 (at A, 
{ x) 7 
re wake (16a): 
Ne pean in which 
2, Py : : : 
eo | Po) = [b2 + As(x-%,)] 4 
tu 
3M 2 2 
7 Tis A ra is ex Pee 
1, | 0 X es X 
Am Ng t J 


ree a 
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2 
c, fat), aR Be be ees -(3 7s 
1 Vain f 7 RT, | aP jo Fo 
a : r=Ry 
: ooX | 
{™,\ ; eR HL G +1) : _ an, ry 4% ZR PN 
a Noi yn MH, | 6 ; 22) ’ 
r 
: P, 7 | <k 
" “ 7 P- ms j ‘ (16b) 
se (: . Mosse ook { in which 
“oop ) 7 a UP 
tt, ro ot ages : K 
P } p2 ] r 
Nehe P(e) = {Fi + Ba be 
a 
(2% W, 
“ -(#) Ty | PP | 
rR, B, nies R | 
' . R, J 


Solutions of equation (13) must satisfy equations (14) and (15) as well as the boundary cordition 


p(X, t) = 0 (17a) 
[ p(R, t) - 0 (176) 
The first «pace derivaiive in equation (13) cas be elimiaated by means of the following 
substitution: 
yo 
pe > = — (13a) 
P, z* 2 l 
fs ates (18h) 
rPy a 


Introducing the new variable v into equations (13), (14) and (15) we take the Laplace tranaform 
under the assumption of initia! rest conditions. 


10 


AE 


THE FRANKLIN INSTITUTE © Laboratories for Research and Developme it 


T-A2049-12 
‘Vhen with 
VCs) v : { se at 
o 
* th) she iat dt 
= fx : 
si(f) = fs fe~™"' dt 
equation (12) transforms ts: 
(19a) 
(19b) 
in which 
12p 1a {x-x, 
eo 3 9 e ~ > “2... 2 2 F (20a) 
He Ao Ue \Pro 7 Pat 
2p : 
1s 4 (20b) 
H, 
F.quation (14) transforma to: } 
ee OMe Ka? Ae ae | 
me aE (Wyo ¢ f Vey Cte) 
ay oO fu ro Hoi =P) ad z 
i 
Ye 
» bh RY) Ps N sy z 
ms? pa a = Moe f ie dr = f(s) (Ub) 
i, HL, te Jn i | 
YW 
i 
ae - a 
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Fowation (15) transferms tee 
( < 
; 7 {¥ ) eared 9 
Heke) ee attat Cp i e  2Ee 5 (t 7 (22a) 
| Uh, - He Pes i, dz lp? 
! : . * “ k : 
= 4 Vv) 4 =~. 
lee 7 Gaya + [sC,-C,+¢, Be} Red SU can Cs (=) ~ § fuuby 
i H, 5 3 "RB, 4a BoHo or 
! : il, R, Pro Rr 
\ 


wil Deve been replaced by means cf the identity 


Fe 
2P,5 
Ke = BO gee (23a) 
I Ay > 
2 
2k, P. 
C6 - —" ¢, (23b) 
R, 
‘Finally. the boundary condition (7) tranatonnn to 
la : 
| Tipe .e) = 0 (24a) 
V(R.s) = 0 » 3 (24b) 


Furthermore, by integrating equation (1S) between the appropriate limits, we can replace 
the integral in eanation (21). Kar this purpose we have: 


Me v bad (a Se \ 
ape Re pay teh (3) = at) ; (25a) 

a ae 5 H, €S \ de p 2 OF v 2 

bees ro at 


(95b) 


C 
=] > 
2 
=i 
wz 
> 
t 
t 
aloi 
| 
Bye 
= 
a 
lypoy 
2; 
o>“ 
os 
i 
3|= 
‘ 
tf > 
P] = 
= 
| eee | 


. - A sae Ct 
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in which 


ak 


| Pini dr , , . (26b) ; a . . 


The fer a its intestal is expressed expliculy in Appendix b. 


Sportroatar selutian f equation (iS) is: 


. : (274): 


a _ h Haye om 
Yor age ryan (27b) 


If v) and v,, ave iadependent solutions which satisfy the homogeneous part of the differen- eee 
thar equation “£9, shen . : ‘ 


VAATED = ¥ fod a) 
CUZ ec oe a (Ba) 
e. ip? 2) oy (Pty. (7? 
H Pad vey Pa vy Ty yg lee? 
' 


v(M yi tp) - VP) 5, (OD 


pee AA (28b) 


Ole, a bg el oraaeay (i Cae 
vp Rav AD - yp Re RE 


\ 


are also solutions having the following properties: 


O(z, 2) = 0 


re Peg 
OU? Bad = QRZ BZ) « I (295: 
A 


== (2,2) + constantt 
ae 


pe See a 8 a erin 


a 


' The Wronskiau of secand-ardar difforwmtio’ sistivine he enug ce iorm y * fixby +> O ine constant. 
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Xx , 
fe pt 0 | ae 
(oe a eee ee ea aS: . 
a (ro5) constant! } ; 7 


‘ 


The solution af (10) vat satisfies the boundary condition (2-4) cau tien be written as: 


* 
ee as rs é i 
: yoy Pp tagp?e aie 30a} 
f ACME 2) i E OCP, 2) :| (30a) 7 
7 
v T N h "any a ' ap : Cy 
V- AGOR ©). ie RP Od, =r Pr) (20b) 
\ n 
in hich V(s) and W(s) are actermined fron equations (21) and (22), 
V sing equation (25) aud mubstituting into equation (22) the boundary values of \ and its 
derivatives obtained by differentiation of equation (30), we have at our disposal two 
rimultaneour equations for the determination of A(s) and h(s). The “‘two by two’” de- 
terninant farmed from the coefficients of the humogencoas parts of these equations is then 
set equal to zero, The location in the complex plane of roots of ‘this characteristic equation ‘ 
determines the stability of the bearing. The characteristic determinant referred to in the i - 
foreecing is as follows: \ 
Panett et ce Sac ec i 
é 


t The Weonskian of second-order differential equations having the form y” + f(x}v » 0 isa constant. 
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(tk) 0 - 


i) : 
dye dys de , 
Aaa T* EN) 56 yn 


be «4x 


(5) 


y= (HY) de y 
oe : 
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2. Stability boundary — Evaluation of Critical (Limiting) Parameters 
Further progress in the analysis depends on a knowledge of the sclutions y, and y, 
the homogeneous part of ¢-yuation (19). Fven with such knowledge the location of roots in the 
complex plane of the characteristic equation (31) might be very involved. Our interest lies in 
assessing the Limiting values of certain bearing parameters which influence stability. Within 
mits of applicability of our equetions, we are guided in the expected range of their magnitudes: ae 


by results obtained by acans of a-simplified analysis summarized in Appendix 2. Different ap- 


proaches were employed for the rectangular and circular geometries aid ticse are described in 
what follows: 


a. Rectangular Geometry 


The homogeneous part of equation (19a) is 


yvn-esz4y = 0 (32a) 


Two independent solutions, which involve Hessel functions of order +23, can he written as [60]: 


1. a 
yor2e ia zig zs) (33a) 
II ‘ : 


in which + 2 3 curcesponds to I and - 2°3 to Il. 


a : 1 $ 
The series representing the Besse! function J,(2iA 2) can be written as [61]: 


HAdin2) o Ga Dy | : a 
AQiA% iA7)y } ——— + ———_. + -—__ oe 
; Meal) * Tied) “Srtysay * ve! 
* With 
1, 3 | 
ve +2 and 2iA” = 3 2/4 ies | 
the first four terms of the expansion of 1 ‘and Yy; become | 
2% 
| (3 ives) | 
kee 3,,2,% Lj 42, 2 ,"2,,,58 
ys T(E) 2 Z+ 5 Z2 “(ées) + 30 Zz (es) + 1485 z (es) (35a) | 
3 i 
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—% 
y “3 : ; é 
(3 7) a ee ee 
Reha resis 1a fe Gai 4 22 tay 2? ay] (36a) 
it ( i’ 3 9 * ane » 367 - ; ; 
r(4y : pear | 
Vos 


‘These truncated series are then substituted into the expression Q(z, ) in (25a) and the 
derivative (dQ a4é)(z, ¢) evaluated with z and ¢ assuming the boundary values of P 2 and 
Bo in the order indicated in the components of the characteristic determinant (31a). Ja forming 
the expressions for the afore-mentioned derivatives of Q, we retain in all subsequent multiplica- 
tions, divisious. and so on, terms of order s* and lower. The assumption that higher pewers in s 
. may be neglected, becuuse of the relative smallness of their coefficients, is justified post factuu 
by noting the: magnitude of the complex-conjugate roots = + iw of the characteristic cvhic to 
which the determinant (31a) is thus reduced. The coefficients of this cubic and those obtained 
by means of the simplified analysis in Appendix 2 are not the same. The effect of this difference 
in results is best illustrated by reference to Figure 3 and the table of results in Appendix 3. 


h. Circular Geometry. 


In this case computer soiutions for a set of discrete values of 3 = iw on the axis of the 
imaginarics were obtained. The probable range of values for o: was dictated by predictions 
based on the simplified analysis and by taking into account the difference in the results of the 
tw analyses for the rectangular case. Stability was assessed from physical conditicas rather than 
proved. If we write 


: 2 
1 ~ (5) 
7 r i2peoR 2 4 . p . P, 
€ = me ee SB is BP - 
2 ‘ 
Kh, Pou ? log, 
ir 
then the homogeneous pirt of equation (19b) becomes 
| : ia aig 
yo+|—-—————]y - 0 ' (37b) . | 
| 1? Tavera | 
t 
che solutions of which | 
De a : i 
y, = 6, (E+ id (OG ; ~ (38b) | 
LH I il l 
! 
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satiefv the conveniently chosen boundary conditions 


y,1} = 1; yy () = 0 


aw 
eo age ee edge 
(1~R log, 3" 


faxes) 


(laplag 2)" | 


and z:he boundary conditions (36b) become 


$()=1; $() = 0 


Sy) = 0; oy 


ul) = 4 (1) = WC) = HL) - 0 


(1) = 1 
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(39b) 


(40b) 


(41b) 


Solutions cf equations (40b), satisfying the conditions (41b) were obtained for a fixed value of BR 


and a set of values of the parameter 


by means of a digital computer. Those shown in Figure A7-1, Appendix 2, correspond to the 


value of a = 0.093, ~ 


12paR? 


+) 2 
Polls 
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eae gow ine position to express numerically the derivatives 
et ob UR) eb AR) 
ae qtr gr ] oh Ee aii . 
Ss ' oy fl) = bey HD 
; of iR) +i det (RD 
“inn. — (42b) 
apie ot Gaye (He iy (RD 
)} 
a aepiec Stee 
dy eg (RD t Ve) CRD 


fer discrete values of cr. 


Pet the bearing eoometey and tie parameters PP. Wh. y and be specified. One 


may tien express the aup width. I. as a function of the supply pressure, P_, by simply equating 
the outflow throush the slot and the inflow through the nozzle. When explicit expressions are 
substituted for quantities which constitute the elements of the determinaat (31b) and a is also 
specified, the only undetermined parameters are the recess depth, 5, and the sugply pressure, P,. 


li fact. the chasacteristic equation assumes ite form: 


afp-o(e.P 1+ ib[s-e,(a. Pi] = 0 (4b) 
inwhich a and b are numerical constants . ‘Therefore, given a@, the problem reduces to finding 
specific values of Po and 6. for which the real and imaginary parts of equation (43b) vanish 


simnltaneuust. e.g. 


&- 2, 7 Ro = 8 (44b) 


ct 


The experimental counterpart of this procedure is to maintain the ambient and recess 
pressures constant and to set the gap width by varying the supply pressure. ‘The recess depth 
is then increased and it» value on the threshold of instability is noted. 


Critical recess depths und supp.y pressures were determined by trial and error. The 
necthod is illusteated by a numerical example in Appendix 4, and in the correspending Figure A4—1 
for a = 0.050 4 plot of g, and ¢, against P,. showing the points of intersection of loci 


curresponding to a = 0.093, is contained in Figure 4. Results of the present analysis are shown 
and compared with those based on the lumped-parameter theory in Figure 5. 
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A Macusston af tresetical Mesults 


Hest of Waving several important parameters on the stability of tae rectangular 
hewn in btxore t Here one finds reasonably good agreement between results based 
cbadted aad bus ped-parameter theories. althoueh the latter theory tends to predict 


vars of uluwaite recess depth which are cither conservative. when the pressure drop across 
visrelatively stall, op slichtly optimistic. when the pressure difference is 
tefeased. Tie aueeeiwert os particularly good for values of the Pressure parameier betucen 0.10 


web thatis far pressure differences in the erder of 19 ped. 


His sere that when the mcess width is decreased by a factor of wo, the allowable 
recess cepth is cearh doubled. so that the critical recess volume remains essentially unaffected 
tna ES perceut increase in stot length. (Conmpate curves Land H in fisure 3.) Values of allow- 
able cocess depih. approvunately tice as large again, may be obtained if the supply nozzle 
diameter ts 


increased by a factor of 3. Compare curves [and JH in Figure 3.) 


This appropriate te remark at this point that decreasing the recess area and increasing 
the surnts nossle diameter have the adverse effect of lowering the load capacity and static stiff- 
hess of the bearing, so that stabilization is often a compromise in the selection of impottant 


design parameters, 


\ aumber of calculations were alse carried out in order te investigate the effect of bearing 
rss on stability. An inersase in mass by a factor of ten showed that the effect was generally 
ceglicible for those rectangular geometries which were, considered. The effect was greatest for 
the configuration corresponding to curve II. Figure 3. characterized by a relatively long slot and 
utd a relatively small supply nozzle. Ata value of the Pressure parameter of approximately 0.04, 
for exemple. the tenfold increase in mass had the effect of reducing the allowable recess depth 


ae 
be 05 percent, 


On the other hand. as shown in Figure 5, the mass had.a significant effect in the case of 
the circular bearing geometry, particularly at low values of the pressure parameter. \t!}) —P / 


ia) 


,, © 0.14, for exomple. an increase in mass by a factor of 3.25 renulted in a 50 percent decreae 
in the allowable recess depth. ‘The deviation cf results based on lumped-parameter theories was 
also significantly greater than for the rectangular configurations, particularly in the case of the 
smallest mass. (Curve Hil. Figure 5.) 


The stualler the mass of a bearing. the greater the rate of response of the bearing tends 
to became. When the response rate is large. the gas in each volume clement hecames partly 
trapped. rather than squeezed out. This. in turn, causes the film “‘stiffuess” in the annulus to 
increase and permit greater “compliance” of a deeper recess without compromising stability. 
When the mass is large, the changes in presaure and density approach the quasi-static condition, 
and the foregoing effects are reversed. These two cases are illustrated diagrammatically and in 
a0 evaugerated manner in Vigure 6 for the conditions of a heavily (a) and lightly (b) loaded 
bearing. In the circular bearing the surface corresponding to the annulus constitutes 96 percent 
of the total aces and the Utst vain eigen luwatd the recess periphery where the gis denaity is 
highest. [f the boundary motion is relatively rapid, the gas cannot readily escape from the slot. 
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Eecanse of this. the infiuer+« of mass on stability is more pronounced in the case of the circular 
bearing. as mat be observed be comparing curves 1 Mand Hl in Figure 5. It can also be seen, 
thit. as the mass increases. the agreement betveen results of the present and simplified analyses 
it preaves, because changes occur nearly quasi-siatically. 


‘The <umplitied. inmped-parameter analysis is based on assumptions af quasi-static chanues 
in pressure and vas cortent and does aut properly account for the effect of masa an the rate of 
tesnouse. fa fact. when the inequality (\2-7). Appendix 2. is written explicitly: 


Ay st As A A AV, Ay ae 
a a Ay : m A, 


the onantity m cancels in the deaominators on both sides, so that the stability condition. is in- 
dependent of mass. Phough perhaps deficient from the quantitative point of view, the simplified, 


Lniped-parometer asalvsis. \ppendix 2. does give explicit relationships between other parameters 
which offer’ 


ieilows: 


a) ? 
and, ay 82 7 Bale 


a ’ (46) 
ll — iW, = Vy ho 
oP, alee = 
"oO 
in which fora given bearing geometry, 
fC > pop 
ee ere ee eae 


My): e(T. PP POD 


It follows that the stability condition also can be written as: 


3) ) 
di dl 
ets S, io (47) 
di dil Pe ee . 
weit, pry he 


The physical interpretation of the above inequality is as follows. ‘Vhe left side of in- 
equaisty (17) represents a change in recess pressure caused by a rapid change in gap width. each 
volume element acting as a closed cylinder in which the gas becomes trapped. ‘The right side of 
inequality (17) represents a change in pressure caused by slow, quasi-static displacement. the 
etuutles Living o flow regulating, rather tian gas-trapping effect. 


vy 


the stability of the bearing, The inequality (42-8), for example. can be written“as” ~ 


ee ee ee ee a a 
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4. Conclusions and Summary of 7 heoretical Pesults 


Nn analysis was nade which teeats the flow within the gas film in a thrust bearing on a 
contiows s. rather than “lumped.” basis Ut shows that within the limitation imposed upon the 
desiza of the bearing Iv other considerations, the following paremeters or combinations of param- 

ters should be mimmized or mivimized in order to insure stability: 


aes 


Winioised Vaxintized 
Derth uf pickets 3 1 Sunply nvzzle diameter 
Siffereace between supply and recess >. Tenvth of annulus 
ues , 3. Area ratio of annulus and pocket 
tb fective mass of bearing regions 


Cne nav also conclude from results of earlier, simplified analyses (2,3.6] which show 


good qualitative azroement with the present one, that high gas temperatures and adiabatic condi- 
lions faver stability. 


‘The most important factor influencing the stability of externaliy pressurized bearings is 
the interaction of the total gas-storage capacitance with the flow-regulating characteristics of 2 
fixed supnly restrictor and a variable annulas. The ratio T. mg. where [. is the load and ~ the 
mass of the bearing, also influences stability. Deeper pockets can generally only be tolerated fre 
very moderate pressure drops across the supply nozzic. ‘The limiting depth increases with she 
nov7le diameter and decreases with the length of the ainulus. as well as with the mass of the 
bearing. ‘These trends are ilastrated in Figures 3 and 5, 


The advantage of recessed packets, channels. ur grooves is that they increase the load 
capacity and static stiffness. The recess perimeter also provides a larger area of admission into 
the annulus. thereby reducing entrance velocities and preventing shocks [19.20.21]. These ad- 
vantages. however. must be weighed very carefully in the light of requirements dictated by 


stability conditions. 
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Wo RNPERIMENTAL APPASATUS 


1. By aren of baperimestad Spparaces 


The main objectives in desi ning the test apparatus were the following: 


a) fo insdre maximum flatness and rigidi:, as well as hardness and material stability, of 
the thrust plates, 
ht fo provide mears of precise alignment. and of accurate methods for checking the par- 


allelism and separation of the bearing surfaces, 


OO) To recteaia the sede of os ciation to an axial direction with minimum friction. 


d) Yo place a yariable+volume, leak-proof cavity ia the center of one thrust plate and to 
‘provide for precise adjustmert and accurate measurement of its depth. 


e) ‘To tumish means of applying a cons‘cnt, non-massive, axial load which would not be 
affected by the motien‘of the bearing. 


muescy band width of oscillations, 


«) To provide sensitive vibration pick-ups and associated electronic and photographic ap- 
paratus in order to detect, measure and record amplitudes of oscillations which are 
small in comparison with the magnitude of the equilibrium gup width, 


h) To furnish a well eeculated, relatively clean and moisture-free air supply. To select | 
icliahle pressure gauges, minometers and flow meters for the purpose of determining 
the static characteristics of the test bearing, 


The tent bearing assembly drawing is shown in Figure 7. The sturdy construction of indi- 
vidual parts and of the entire structure, particularly the plates and columns, ie self-evident. All 
members were machined to extremely close tolerances, subjected to proper heat treatment and in- 
spected prior to, during. and after manvfacture and assembly, The upper thrust plate On cculd be 
aligned | me ins of three 2 inch diameter, 48 threads per inch, adjusting screws (2) and securely 
bolted to the supporting columns & - The lower thrust plate was integral with the shaft ) : 
Hoth thrust faces were chrome plated and lapped flat to the erder of 1-1/2 interference bands of 
monochromuaiic sodium light (approx. 20 x 10°6 inches). Final inspection, using 6 inch diaineter 
opticat flats, showed that the convexity of a symmetric and almos? circular central region of both 
plates did not exceed one interference band up io approximately four-fifths of the 5 inch diameters. 
Three sensing probes @ot an air gauging system, incorporated in the upper plate, were equally 
spaced along the perimeter of the thrust face und lapped together with it. 


An exact determinution of deflections and modes of vibrations of the thrust plates was not 
attempted hecause of the complexity of the problem involving thick plates, supported in a manner 
described in the foregoing. A conservative estimate, based on approximate plate theory, the 
a a Ne TE 
t Cireled numbers correspond to those in Figure 7 and refer to assembly componenis, 
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. ; 
iiple support aiong the column pitch circle diameter and a concentrated 500 !b, 
s (foe bend Slowed this the dedlecrion of the center of the upper thrust piete with reference to its edge 
© RE See ae he idee at 20 x tu inches. The deflection of the lower thrust plate. in which the canti- 
ime riag was subjected to pressures which in that region approached the ambient 
velbe. would be negligible in comparison. It was concluded, therefore, that load deflections 
swonld have a Compensating effect on the slight, initial convexity of the surfaces, 
: fe a ' 


“assumption o 


tecrptibe ls apverty 


«The theary of vibrations of plates, even in apnroximate form, is verv complex. On the as- 
suf ation that the frequency , o,. of the lowest mode of vibration could be estimated from the 


“ermyrtion Fae t: 


| ( Ee 7 
ade ar ag ee) 


; ‘in which -a@ is a constant depending on the method of support, Eis Young’s modulus, » the 
; Poisson's rating p the mass density of the metal, r and t the radius and thickness of the plate. 
“ahe fowest frequency was estimated to be of the order of 104 radians per second and, consequently,. 
céasilecably higher than the anticipated frequency band width of auto-oscillations. 


~The upper thrust plate, Figure 8, contained a centrally located recess, the depth of which 
could be varied by means of a sliding piston © and a tiicsometer adjusting screw @. Three 
tadtly. fitted “O" rings served us sliding seals. The piston cauld be locked in any desired posi- 
tion and the recess depth varied from zero te 3/4 inches. Air was supplied through a 7:16 inch 
| diametce, drilled duct in the upper plate. It flowed into a chamber, the volume of which comprised 
Z | the seated space between the ‘*0" rings and the bore of the piston, Figure 9. The flow restric- 
tion between this chamber and the recess was effected by means of a short, bell-mouthed, 0.055 
\ inch diameter nozzle @).t The drawing pf an auxiliary piston (Figure 10}, shows a 0.0135 inch 
diameter pressure tap, located ai i 1-32 {aches from the piston center, The tap was connected 
through an internal duct to a gauge which indicated the static, equilibrium recess pressure at that 
point.. In the modified version of the piston, the pressure tap was centrally located and plugged 
when dynamic tests were conducted, 


‘The drawing of the shaft and lower thrust plate is shown in Figure 13. The shaft was hol- 
low svthat the movable mass corresponded to 2% lbs. The weight could be increased to 45 Ibs. 
by attaching a 20 1b. bronze diac to the shaft, midway between the 3 inch journals, (It will be re- 
called that theoretical stability loci were constructed for masnes corresponding to weights of ap- 

., proximately 13.9, 25 and 45 bs.) Rigidity requirements,-however, precluded a less massive shafe 
“and thrust plate constuction and no provision could be made to obtain experimental data for the 
smallest mans. the «inch diametes, ¢ rcular area of the plugged end of the shaft provided a sur- 
face exposed to the pressure of the loading cylinder 9 Ly means of which the thrust load was 
applied to the bearing. Impulse turbine buckets were milled in the shaft to provide a high speed 
drive for experimental work invelving rotation which could be conducted upon the completion of 


a er ee 
VIn a modified, later version of the piston, provision was made for replacement of a Single racteiotar hy 4 Oo 
ori2 nozzles, In eech case the total crossesectional throat atea was adjusted to give approximately cal- 

culated rates of inass flow (see Figures 11 and 12). 


[4 
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this ph ise of the .avestigation.t The size, concentricity, roundness and taper of the shaft jour- 
tials were found to fal! within tolerances specified on drawings. The total indicated runout of the: 
trust face with respect to the journals, mecsured on a 5-3/4 inch diameter, was less than 0.0002 
inches. The shaft was bataaced to within 0.00095 inch-pounds. 


‘The problem of restraining the motion of the shaft and the lower thrust plate to a purely 
axial direction by means of low-friction guides was solved by mounting the shaft in extern ally- 
pressurized, air-lubricated journal hearings @3 . There were two main reasons tor this tvpe of 
later} sunnort, First, the analysis was predicted on the assumption that external damping was 
nestivihle, since primary interest was only in those stability parameters associated with th: hear- 
ing itself, Secondly, from the practical point of view. it was desirable to simulate as far as pos- 
sible a method of gas bearing support which could be eventually applied to small turbines, com- 
pressors or centrifugers. A drawing of one bearing is shown in Figure 14. The main feature of the 
design was a set of six, grooved, centrally located and equally spaced pressure pads. Air was 
supplied to the pads through 0.016 inch diameter, brass orifices () , shown ia Figure 15. The. 
purpose of the grooves was to minimize the storage volume of each ad, while extending a rela- 
tively uniform, high-pressure aren to the outer perimeter of the grooves. The net effect of this ar- 
rangement contributed to the stability of the bearing and produced a steep joad-dellection charac- 
teristic, Figure 16. Needle valves were placed in drilled air ducts leading to alternate pads, 
The function of the needle valves was t allow for controlled misalignment of the shaft within the aes 
limits of the bearing clearance. By varving the pressure in individual pads, am excellent inethod ae 
of making final adjustments in the alignment of the thrust surfaces was provided, furnishing im- oe 
provement on the high degree of parallelism already achieved bv means of the levelling scrcws: 


y 
i 
y 


The bearings were dowel-ninned and bolted to the supporting members pricr to final ma- : : > ie 
chining, during which both bearings and the sleeve of the loading cylinder were yround to sizé in i 
one operation. The grinding operation was performed with a degree of precision which made it 
possible tu span all three bores with a single plug gage. barely 0.000% inches smaller than the 
niniaium, measured bearing diameter. All grooves were indicated or the grinding machine and Fs 
their depth varied from 0.006 to 0.007 inches. The orifice inserts, needle valves and the »noular ' gat 
supply chamber of the bearing were sealed by means of ‘‘Q”’ rings. : 


It wee estimated that a bearing stiffnesntt of 3 x 105 to 5 « 105 pounds per inch could be 

realized when the journal was in the concentric position. With the hearings and the shaft assumed 

as a spring-mase system loving two degrees of freedom, approximate values of frequencies of j 

natural modes of vibrations were calculated by considering a pivoting motion of the shaft about ’ 

one bearing, the other journal being acted upon by a spring. The “‘constrained” frequencies, p, 
| 
} 
i 


ate given by (58): 
k 1/2 
pet fe oe 


tThe grooves at the end of the shaft were intended to suppress inntability (whirl) due to hydrodynamic 
forces jn the clearance of the loading cylinder slerve. , j 


Tf'the stope of the load-deflection curve is sometimes referred to as the “‘atiffnesn’’ of muse correctly, an i 
the “‘static atiffnenn’’. (See Figare 16.) 
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L. the distance from the chosen pivot to the 
center of gravity of the shaft, k the stiffuess, m and [. the mass and moment of inertia uf the 
shaft about itz penter of gravity. ‘The estimated natural frequencies far the 25 ib. shaft were in 
the order of 3.x 10 and 3 x 303 eadians per second and could be substantially reduced, or raised, 
by ving the supply pressure, [t was concluded, therefore, that in the absence of appreciable 
unbalaeee and movents in planes passing thrusgh the shaft axis, the pitching motion of the shaft 
within the clearance of the journal bearings would be relatively small ia comparison with the am- 
pinude ot axial oscillations. 


inwhich « is the span length between the bearings, 4 


The loading evlinder Ca consisted of a 7-1 2.x 3-1/2 inches pressure vessel attached to 
the lower end of the supporting steuctuce. A constant pressure maintained in the cylinder and act- 
ing on the ead of the shaft protruding throug the evlinder slcave, balanced the load applied to the 
lewer thrust plate. Fur snail amplitudes of oscillation, volume changes would be negligible in 
comparison with 1! hic ine hes of cylider displacement. The cylinder pressure would remain 
essentially constant, unless the frequency of oscillation produced a condition of accoustic reso- 
nance in the chamber, in which case the shape and volume of the cylinder could be altered accord- 


ingly. The evlinder provided « simple method of applying a coustant, non-massive load to the test 
bearina, 


the lower ends of the 3 incn diameter ¢ turns of the test rig wese prens-fitted into a 2 
inch thick plate (9) in which one of the jowma! bustin gs was iaousied. The upper ends of the 


oe oory roan . 
colunns were located and bolted to a similar plare G5) ecsleining the second journs] bearing. 
~ 


The column locating holes in both plates were jig bored in tandem. The cuuccairicity and square- 
ness of the columns and locatiug bosses were indicated from the hoad of iio precision grinder aft- 
er the bearings and the bore of the loading cylinder slecve wore finished to size. ‘The assembly 
was bolted to a.massive and very rigid tripod, The weight of the entire structure was approximate. 
ly 750 lbs. Inthe course of experimentation it was found necessary to interpose relatively soft 
isolators between the tripod and the foundation floor, The reason for dsing so is discussed later 
in the text, A schematic of the air supply system is shown in Figure U7. Photugrapha of tear ap. 
paratus, its components aad of auxiliary equipment are shown in Figures 18 through 25. 


2 Instraments -- Pefursance and Calibration 


a. Air Gauges 


The threefold function of the air gauging svstem. briefly referred to in the foregoing sec- 
tion, was (al to facilitate the alignment ef thrust surtaces to a high degree of parallelism, (b) to 
measure the equilibrium gap width, Ho, and {.) to provide a quick and direct method of calibrating 
the primary capacitance probe without disturbing its initial position aad offset. 


The well-known method of measuring small displacements in terms of pressure utiiizes two 
fleny resiricters, placed in series, a constant, well regulated pressure source and suitable pressure 
indicators, Air eaters an interme:tiate chambe- from a constant pressure reservoir through a small, 
fixed (or adjustable) orifice. ‘The variation of pressure in the intermediate chamber is recorded 
and depends on the supply pressure fevel, as well as the relative size and flow characteristics of 
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the cesteictors, The cross-sectional flow area of the second restrictor is made to vary -vith dis- 
placement, so that resistance to flow is fumished by a slot, ur a circumferential orifice formed 
between the surfaces of the sending probe and the displaced member. 


bora fixed orifice size. the slope and curvature of the pressue-displacemert ~haracteris- 
tic depend mainly on the size and geometry of the prohe restriction, and on the magnitude of the 
suppis pressure. ‘The characteristics are nonlinear and exhibit an inflection point, in the neigh- 
horhosd of which. the amplification is highest. The probe can be made in the shape of a short, 
thinswafled or a thick-walled cylinder. In the lattcr case, the restictioa between the probe and 
the displaced member tus che form of a narrow, redial slot rather than tnat of a circumferential 
orifice, Alt ether factors being equal a slot-like, variable restrictor produces a higher sensitiv. 
ity. thought ai some expense of Linearity an the pie ssure-displacement characteristict, This effect 


is due toa more pronounced infisence of viseaus torces in a natcow slot than in a circumferential 
orifice, : 


It was considered that a low supply-pessure gauging sv stem, inte gral with the upper 
theast plate, would have an advantage cover commercially available units, wich generally operate 
at supolv-precsure levels of thirty to sixty pounds per square inch. The system was ts be e-aven- 
contly incorporated into the test cig assembly, with probes rigidliy mounted in the upper thrust 
plate and lapped in the plane of the thrust surface. A well regulated supply head of two to three 
feet of water could be obtained by means of a standpipe, with the excess air in the supply reser- 
voir being allowed to escape at the end of an immersed tube. [t was also desirable to avoid the 
inipingeaunt of celatively high velocity probe jets associated with high pressure systems. since 
this impingement might hive involved difficeter oy ty valibrativu, aod possibly minor effects of heat 

‘transfor on the accuracy of measurement. 


An estimate of orifice and pishs dimensions for a low pressure sysiem was made by assum- 
ing the flow to be incompressible. Dimensionless pressure-dieplacement charecteristics were ob- 


‘wiered either as (a) a narrow, radial slot or (b) a circumferential orifice. Results af these pre- 
liminary calculations used in establishing approximate values of design parameters are shown in 
Figure 26, which also contains a schematic of the supply head regulator. © 


The effect of changes in barometric Pressure on gap width measurement could be estimated 
from equations (8) and (9), the flow through the orifice and probe restrictions heing esseitially 
analogous to that through the bearing supply nozzle and gap. The gap width is then given by- 


(P,P 


Ss ty 


Pe ~ es t 


~P2yV? 


W3 = K 
° 


a 
t The princiolns ts 


ener enanestncneneniniiennstaees 
r rodin Jd taniuiug the prensuredispiacement characteristic of an air nauge are esnen- 
tially similar ie the-e which apply io le xd-displecement characteristica of extemally pressurized, au 
hearings. 
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immueit Fy, Pe and BP, nov reler to the absolute values of the supply, intermediate and am- 
haent pressures of the air gauce system, and K isa constant. Taking the logarithm of both sides 
of thus oxpression and differes tiating, the aypruximate error in gap width measurement, associated 
with a vaange of t VP oan barecmetric pressure, ia: 


Ally |. Att 


[~_— 


Ho ~ P 


where inthe last step it was assumed that the pressure differences are sensibly small, e.g., 
ot DP ~ fee > P.A change of + 0.5 psi, relative to the barometer reading during cali- 
bration, would cause an error of JF 0.5 percent in the measured gap width, less than 10 n-icro- 
invhes, 


Initial difficulties were encountered because of supply pressure fluctuations caused by the 
surface tension of bubbles, forming and breaking away from drillings at the immersed end of the 
tube. ‘This was circumvented bv inserting a damper in the form of a series of perforated discs in 
the reservoir end of the tnbe and reduced the ripple to an almost imperceptible magnitude, relative 
ta the smallest manometer scale reading ot 1/10 of an inch. The supply pressure was adjusted to 
a steady, constant value of 39.2 inches of water by varying the setting of a fine, long taper needle 
valve, 


The calibration wan performed in steps of 0.0001 inches over a gap width range from 0.001 
to 0.002 inches. Square, 1 x Linch precision gauge blocks, having a dimensional and parallelism 
tolerance of 1 microinches, were used for this purpose. Two blocks, differing in thickness by a 
multiple of 0.0001 iaches, and a third block were wrung together. The step gauye created in this 
fashion was then wrung tugether with the lapped thrust surface, the step extending from the edge 
af and overlapping the entire probe diameter. (See schematic in Figure 27.) The probes were 
spaced 120 degrees apart and located 1/4 of an inch from the edge of the upper thrust surface. 
The blocks were taken apart and reassembled after each round of manometer readings for all three 
probes, This was repeated 21 times for each probe and everyone of the 11 calibration points, a 
total of G91 manometer readings. The gauge blocks were assembled at random and it was reason- 
able to assume the arithmetic average of 21 manometer readings as the most probable value cor- 
responding to the nominal dimension of the step gauge offset. The curvature of the pressure-dis- 
placement characteristic in the experimental range of gap width was very slight. The maximum 
difference in any set of 21 manometer readings was less than 0.4 inches of water. Ata represent- 
ative sensitivity of 0.15 inches of water per 10°5 inches displacement, this was equivalent to 
2.7 x 10°5 inches, or. Jess than 3 percent of minimum gap width in the experimental range. Lenst- 
square-fitted calibration curves, Fivures 27, 28 and 29, were considered adequate for gap width 
measurement to an accuracy of 1 to 2 percent. 


b. Vibrometers = Primary and Auxiliary Probe Systems 


The validity of assumptions made in the analysis increases if the amplitude of motion is 
small in comparison with the equilibrium gap width. For reasons arising from practical considera- 
tions, such as the impossibility of insuring perfect isolation, ideal rigidity of structural members. 
absolute steadiness of air flow etc., the feasibility of detection and, above all, meaningful inter- 
pretation ct amplitudes in the order of 10°° inches was duubtful. It was arbitrarily decided, there- 
fore. to regerd amplitudes of oncillation in the order of one percent of the gap width, or to 10°5 
inches, as small, : 


2» 


SS 
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The fupetion af ea h probe svete sas te biter ot dvaamie capacitance changes in ters 


oftelitive disphiceneur ofthe ae dbhe and the vibrating member, Sie primary prod sy stent, stow 
onthe dedt si hoof Pieure SO was qeed in measuring the avial oseithations with which the experi- 
eek tad caw ea ile a een, ola await profe svatem, shewn on the right side of Figure 30 
sraploved bodpactive pickeup, shiek could be mounted in several positions along and 


aeagind the slaft. Gre tection of dre auvilian probe wart. detect transverse males of vibrations 


; wh ri byevstt Le coupled with the ayial motion of the shaft. It wis, of course, desirable to reduce 
‘ ; a s hes thodnotion tea hare qinionmes: ond dhe airelubri cated aide bearings were desizned with this ohjec- 
tive unwind The auvilian sestem thas piasea a moritoring cole in the course of experimentation. 

The aperatim ofthese vihremeters was bised on different principles. Their operation, perform~ 


ance caldbration ae deseribed in the next paragraphs. 


Ie thas wehe oovagrin, ofthe aimen instrament. a Fielden Proximity Meter, model 
TM-OS1-1. the variable probe capacitence. Cl. is shown chanted across one of the two balancing 


we excitation was a 450 K.C. oscillator 


ctpacitars ofan impedance bridge, The source of brid 
(modulated by Sp Ke owhen using the weter for static measurements). The out-of-halance voltage 
efthe buidge was applied in succession te the tube erids of a twoestage, radio frequency amplifier. 
: The gum of the amplifier contd be conteolled by a variable. tapped resistor in the cathode circuit 
zs of the amplifier tubes bx means ofa sixeposition switch. providing an attenuation up ta 50 deci- 
hele in steps oof LO decthbels. tte bridge could be dalanced by means of capacitors, C, ry 


and the potentianeter Ro woich balanced the probe lead resistance. The free 
instrmeent. Fiaure Sho was flat from DLC, te approximately 1200 cps and. therefore, more than, adc- 
Dmeastrements ithe anticipated 100-100 cps band width, The ouiput of the insirument 
wascapplied ty one pair of tecmingds of a dual trace preauplifier of the cathode say oscilloscope. 
The preamplifier had cine sensitivity settings ia the range of 0.05. .m te 20v cm and a variable 
veemeradjustuent, The frequency cesponse of the oscilloscope-preamplifier unit (Tektronix, type 
SQ and dualetrace preemplifier. tepe 53 5404 was Yderibels down at 5 megacycles, considerably 


aiiate f- 


higher than that of the Fielden ‘nstrument. 


The primary jeobe. Figuce 24, consisted of al tinch diameter, bras» disc, mounted in an 
insulating sleeve. The dise any sleeve assembly were pressed inte the opening of a 4.5 inch long 
probe hotder and Lapped square with its axis. The bedy of the huider, which also served as part ? 
ebte the upper thrust plate by means of a split collet and a 


: of the outer shield i : 
: : «gripping clamp, Wich the radia frequency amplifier contra’ adjusted to approximatels one third of 
, eamplifier vertical deflection setting of | volt per centimeter, the sensitiv- 


gain and 
its corresponded to 1 centimeter vertical deflection per 10° inckes displacement when viewed on 


mavin 


the disphas of the cathode vas oscilloscope The observed noise level and drift at this setting, 
after an iaitiel, one hour wari up perivd, was of the order of 0.05 to 0.7 valts. The maximum room 
temperatnre variation during eet tect run was between Thand TT degrees Fakzenheit. The probe 
lead was double-shielded tu prevent stray voltage pick-up and relatively short, in order to mini- 
mize capacitance in parallel with that of the probe. For this reason also, the thickness of the in- . ! 
sulating sleeve (diciectric constant of 2) was U.1 inches and, therefore, large in comparison with i 
the initial, 0.00% rach offset between the | Linch probe disc and the lower thrust plate, i 


Following the pareliel aligaaeat of thrust surfaces, the probe was calibrated in three posi- . 
tian enrtere Lippe cateaiels ine cad and tidpeints of lie experiuental ragge of gap . 
widths, The calibration was performed, using the aie gauges as a standard. The calibration ‘ 
; . 
an : 


me, 


ae et ! : ee 7 && 4 
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curves are shawn in Figures 32, 33 and 34. rom their appearance it could be concluded that for 
amplitudes of disolicement-ia the order of Lx 107 inches to 2 x 1075 inches, no large error would 
be Greurred by assuming linearity in the reighboshood of null points aad a sensitivity value of | 
volt per 1079 inches displacement, 


The schematic diagram of the auxiliary instrament (Decker, Delta Unit 902-1 and Monitor 
Tait 902-2) shews the capacitance probe, C,. mounted in one of the positions along the shaft. 
Tic main component of the system was a tube, containing gast under reduced pressure. The tube 
was placed hetween two external electrodes connected to a 250 KC oscillator. The excitation 
cused a lumiaeus discharge, se that the two internal electrodes of the tube were in contact with 
toe plasma of th» ionized gas. Any asymmetry in the spacing of the electrodes would give rise to 
+ potential difference between the internal electrodes, caused by the migration of electrons. The 
effoct of asvametry could aiso he produced by changing the magnitude of the sensing capacitor, 
C,. relative to the balancing capacitor, Cy, both shown in the schematic diagram of Figure 30. 
The operation of the auxiliary probe, Cy, was based on this principle [62]. A cathode follower 
provided a low anpedance outpet at the preamplifier terminals and a filtering network eliminated 
the fraction of oscillator excitation passed through the transducer. 


The transducer, probe, mounting bracket and calibration fixture are shown in Figure 25. 
The probe and the balancing capacitor were mounted very close to the internal tube electrodes and 
all leads were relatively short in order to minimize the adverse effect of stray capacitances on the 
output and frequency response of the instrument, 


The frequency response of the auxiliary instrument, Figure 31, was flat to approximately 
300 cps. 0.5 decibels down at 490 cps and, therefore, quite adequate in the anticipated, experi- 


mental frequency band width. The calibration curve of the auxiliary probe, Figure 35, displayed a. 


greater degree of nonlinearity than that of the primary probe. The sensitivity in the neighborhood 
of the nuii point was approximately U.UZ voits per 10°" inches displacement, representing 0.4 cen- 
timeters vertical deflection per 10°° inches displacement on the cathode ray oscilloscope display 
at the highest preamplifier sensitivity setting of 0.05 volts per centimeter. The observed noise 
level was of the order of 0.005 volts. Since the maximum sensitivity of the auxiliary probe was 
considerably less than that of the primarv, the monitoring test of lateral shaft vibrations was con- 
ducted under conditions of axial oscillations approximately twenty times as large as the ampli- 
tudes of self-excited vibrations observed in the course of stability experiments. 


- . 


c. Pressure Gauges, tlowmeters and Regulators 


The supply, recess and lvading piston pressures were measured hy means of two, dual 
bourdon tube construction, eight and one half inch dial diameter, absolute pressure gauges. The 
gauges were graduated in steps of 0.1 and 0.2 psi over their respective ranges of O~100 and 


0-150 psia. The ranges were covered in approximately two revolutions of the puinter. The scale 


tad a mirrored surface reflecting the fine blade of the pointer to eliminate parallax errors and the 
size of the smallest subdivisions permitted reliable estimates of readings to the nearest 0.05 psi. 


The neeaaure gaara woes caliheatad hy the manufactuear hy means of 3 precision mercury | 


eter, a duplicate of a Bureau of Standards reference barometer. The accuracy of the pressure 


Tt Argon or Neoo 
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gauges was guaranteed to be within 0.02 percent of full scale deflection, representing 0.027 per 
svat for the 0-100 psia gauge and 0.041 percent for the 0-150 peia gauge at 72.5 psia, the lowest 
test reading. The scales were hand drawn to match over one hundred calibration points. The hys- 
‘eresis, backlash and ‘dead zone’? of the pressure gauges vcve entirely negligible. As an addi- 
timal check, the gauges were conne ‘ted to the sare pressure source and their readings compared. 
No perceptible difference was observed in the range from 50 to 100 psia. 


Flow measurements were performed, using a variable area flow meter (rotameter) in con- 
jun: tion with a calibrated, bellows type gas meter. The position of the rotameter float was regis- 
teved along the scale of the 600 millimeter long tube for any particular setting of the beariag sup-_ 
ply pressure. With this setting maintained and the position of the rotameter float unchanged, the 
flow was diverted to the gas meter. The iotameter thus merely ve-ved as a secondary, indicating 
instrument, while the actual flow rates were determined from the gas meter. readings. The calibra- 
tion of the gas meter was performed by means of an inverted-tank-ty pe gasometer of 20 cubic feet 
capacity, used as a standard by a local utility company, Five calibration runs were made after 
the final adjustment of the counter mechanism, at a representative bearing flow rate of 5 cubic 
feet per minute, reterred to standard atmospheric conditions. (Approximately four revolutions of 
the five cubic foot counter.) The differences in the final positions of the gasometer displacement 
indicator were negligible in comparison with the tota. displacement. ‘The speed of the counter, 
approximately one revolution per mintite, was slow enough to enable the observer to estimate the 
position of the pointer, relative to a predetermined mark, at least within 6 degrces, or better. 
This, in terms of displacement, corresponded to approximately 0.5 percent of the casometer vol- 


ume, 


The «otameter float was operatedin the upper: alfofthe tube. Slight uscillations of the’ float 
were observed in the uppermost tube region, but their peak-to-peak amplitudes were of the order of 
2 to 3 millimeters, as compared with the 600 millimeter length of the tube. An estimate of 2 to 3 
percent accuracy in flow measurements obtained in the course of experimentation appears to be 
reasonable, 


The pressure reculators uscd in the test-hearing air-supply system, Figure 17, were capa- 
ble of maintaining a downstream pressure setting within 0.1 psi over their entire flow range and a 
25 psi variation in delivery pressure. Qne regulator was installed ahead of a 7.5 cubic foot air 


‘storage tank and individual regulators were plaved in the supply lines leading ‘row the tank to the 


loading cylinder, test and guide bearings. The performance of the regulators was most satisfac- 
tory and pressure settings cuuld be accurately maintained, with drift not exceeding 0.05 psi over 
periods of two to three hours. 
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\. EAPERIMENTS 


Loootatie bests 


the calibrated air siuges furnished the cans of bearing surface alignment. The 25 Ib. 
stuli was teserted in the peessurized wide beartags and the upper thrust plate placed on the col- 
wrin becating basses. The conceniriz in of the hearing menhers: wis checked and adjusted to with. 
U.UOT inches. foe thrust surfaces were ten brought inte parallel alignment by means of level- 
line screws and vis sanges and secured tothe coluaos. As an additional check of parallelism, 
the fatawing test was performed. V spacer. aprravimately £,75'x L077 inches thick, was lightly 
clamped between the thrust surfaces at twelee, thirty degree. angular intervals along their perim- 
eter, Manometer readings of sie air zauze only were registered for each position. This was re- 


peated nine times andthe uwithwelie averave of readings for each position was taken. An identical 
test was carried our with a spacer 1.53. 10°" inches thick. In the process of sampling, the spacer 
was efamned bya constant, (13 Ih. force. at an average distance of two inches from the bearing 
center, The guide bearings. spaced nine inches upart, were pressurized and their stiffness in the 
concentric position, Figure 16, was in the order of 1x 10°? pounds per inch per bearing, Any mis- 
alignment, therefore. caused by the moment of the small clamping force, would ke negligible. The 
results of this alignment test are shawn in Figure 36, in which the circularly spaced numbers re- 
fer to average aanoneter readings of wauce No. t obtained with both spacers. .By reference to the 
caltbration curve. in Figure 27. it can be Seen that the thrust surfaces were paratiel within one 
percent of the equilibemmm gap widths, that is inclined under an angle in the oraer of 4x 10°5 de 
grees.? ‘The shaft thrust face had a lotal indicate 1 renout of 0.0002" relative to the journals and 
the alignment would have been destraved, had the shaft Seen allowed to drift ia its guide bearings, 
An additional constraint, therefore, was fumished by means uf two opposed and narrowly spaced 
air jets, impinginy oa a light, thin fin, attached to the shaft collac. (See Figure 18, for example.) 
The alignment was thus performed with the shaft maintained in the same position relative to the 
snide bearings. 


The analysis was based on the existence of an essentially uniform pressure plateau, ex- 


over the entire, one inch diameter recess area. It was postulated that, as long as the re- 


tiseh large, the effects of tie air jet issuing from the 0.055 inch diameter 


nozzle and the acceleration of che fluid ay it entered the slot would be local aad detract little 

from the recess volume conceived as a constant pressure reservoir, Since the stability boundary 
toe . -> i i 

for critical recess depths in the order of 10°° 


tendin 


Teas 


inches was of particular interest, the question which 
presented itself was the degree to which the fore acing suppositions were valid for recess depths 
smaller than the nozzle throat diaineier und comparable in magmtude with the bearing gap width, 
The following possibititiés deserved careful consideration. {a) For &~ 1 4 dy 
crea bounded by the nuzcle diameter and its projection on the lower thrust surfa 
the nozzle throat as the ceatrolling restriction. (h) The recess, its depth sufficiently reduced, 
could act as an inward extension of the bearing gap and virtually cease to function as an interme- 
diate reservoir. (c) ‘The accelerating flow in the diverging region adjoining the nozzle exit could 
attain supersonic velocities, follawed by shock and a transition to the predominantly viscous flow 
a ee eC eee aE 
‘Interterence fringes observed under optical flats indicated thar each 1.5 x 107° to 2.0 x 1975 


convexity waa quite concentric and approaching in shape that of & surface of revolution. 
or inclination, therefore could be thought of in reference to the circular bases of those 


+ the evlindrical 
ce would replace 
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; : Utee spon. ; 
region in the slot f19, 20, 21], a situation analogous to dat which might exist in a duct feu by a ; 


- converging diverging nuzzle, . 


Vtest was conducted, in which the load and a representative value of the supply pressure. 
37.5 pia. were held constant. while the recess was being reduced in depth fram a relatively large 
value. The objective of the test was to determine the “‘cut-off” point, at which any further reduc - 
tion in lepth would cause a change in the steady-state flow pattern and pressure distribution. 
This weald thea be wecompanied by an easily detectable and measurable departure from the initial 
gan setting, Results of this test, presented in Figure 37, indicate that. for a centeatly located, 
Q.055 inch diameter nozzle, the “cut-off? point corcespunded to a recess depth of 6.020 inches. 


dince the critical stability paremeter was the volume. rather than the depth of the recess, 
the fir hipt ata remedy was to surround the nozzle exit by a 6.375 inch diameter, 
6,020 inches deep cavity (see Figure 22), This, however, did not provide an entirely satisfactory 2 
solution, since for small critical volumes the delineation between slot and recess becan.e rather 
saune. For this reason, and because the linear dimensions of ti restrictor influence the degree 
of uncteadiness in pulsating flow, the pistons were replaced hy mailified version, Figure II, which 


a success yl 


incorporated interchangeable, multipte-nozzle plates. ‘The plates, Figures 12 and 23.. contained ‘ 
i either four, nine, or twelve short, bell-mouthed nozzles, equally spaced on a 0.5 inch diameter cir- 


cle, ‘Their total, cross-sectional throat area was later adjusted to give a flow characteristic ap- 
proximating the calculated one within the experimental range of pressure ratios. The twelve noz- 
cle plate was used in the concluding experiments. (ta supply pressure of 47.5 psia, a reduction 
in cecess depth flow 0.008 inches, the ‘cut-off’? point, to 6.006 inches, was accompanied by a de- 
crease in gap width ii the ordec of two percent only. ‘The effects of decrease in recess depth an 
gp width for single and multiple nozzle feeders are compared in Figure 37. 


The calculated load capacity of the test bearing, based on the assur:ptions of laminar flow 
inthe clearance, negligible entrance effects and uniformity of s-essure in the recess, varies al- 4 
most linearly with the recess pressure and is independent of tie supply pressure. 


In order to ascertain the extent to which the actual conditions conformed to this idealiza- 
tron, a test was conducted in which (a) the load was adjusted to maintain a constant recess pres- 
suret of 73.5 psia. or (b) the load was held constant at the calculated value of $72 Tb, and the - 
Tecess pressureT was aliowed to seek its own level. The result of this test is presented in Fia~ 
ure 38 and indicates that in case ‘a) an increase of the supply pressure to 95 psia was acconi- 


panied by a 5 percent reduction in load capacity, whereas in case th) a fb percent increase in re- 
cess pressure was observed. Throughout this test the recess depth was larger than the “cut-off” | 
value, referred to in the preceding paragraphs. 


4 ‘The desirability of minimizing differences between the actual and theoretical, steady-state 
hearing characteristics was of corsiderable importance with reference to the mass flow rate, since , 
its partial derivatives with respect to gap width and pressure, equation (16b), played an important 
tole in establishing the limiting values of stability parameters, Preliminary flow measurements in- 
dicated, that in the experimental range of pressure ratios the overall discharge coefficient, using a 


i 
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34 
‘ KS Tt + 
7 ‘ Lad , 
‘ a i ee ica * B ~ , 
. Le . ‘ *. ‘ 


~~ 


enn ong. 


THF. FRANKLIN INSTITUTE + Laboratories for Research and Development 
T-A2049-12 


iselve nozzle supply restrictor, varied approximately from 0.78 to 87. The total dhroat areca was 


then increased beapproximatels 10 percent, in order to narrow down the gap between actual and 


calculated flow rates. 7 


Mass flow rates were determined in ihe range of supply pressures from 75 te 95 psia, with 
aconstant load of 172 'b., the calculated value corresponding to a recess pressure of 73.5 psia 
and anbieat pressure of LET psia., applied to the test bearing. (A condition analogous to that - 
ander which stability tests were conducted.) The recess depth was always maintained ata value 
areater than the “‘cut-ofi’” point referred to in the preceding paragraphs. ‘The mass flow rates in 
the experimental range of supply pressures, obtained with slightly increased nozzle dianeters, 
were Of to 97 percent of the calculated values, while deviations of measured gap widths were in 
the order of ane percent. Results of mass flow and gap width measurements are presented in Fig- 

aU Mt v 


is 


ures SO and W. The correspouding bearing Reynolds number, Nr * 


it Jaa. 


°o 
wives in Figure 41, Mass flow rates and the Reynolds number obtained in a: supply pressure 
range for a bearing load of 222 !b. and corresponding tc a stability test conducted under reduced 


load conditions are given in Figures $2 and 43, 


It seemed highly unlikely that absolute, one-to-one correspondence between calculated and 
actual, steady-state values of pressure. gap widths, and fluw rates could be simullaneue ly 
achieved throughout the experimental range. It was considered that, in the presence of slight de- 
\ivtions in several static bearing characteristics from those of the idealized model, a fair compro- 
mise was reached, ate! that the deviations were of a reasonable order of magnitude. 


2. Preliminary Stability Tests 

Initial differences in calculated and actual, static bearing characteristics and subsequent 
modifications incorporated in the apparatus, were described in the preceding section. These mud- 
ifications, and a number of additional changes, were also dictated by evidence gathered from re- 
sults of simultaneous! conductcd. preliminary, dynamic tests. The preliminary tests furnished 


section an account of preliminary experiments is presented in a logical sequence, not necessarily 
in the order in which they were conducted. 


The dynamic experiments, with the exception of one comparison test, in which the recess 
pressurett was held constant, were conducted under constant lead conditions. Results of static 
tests indicated that, in the latter case, deviations frow predicted, static bearing characteristics | 
appeared to be less pronounced. Moreover, since the.reccss pressure could be measured at one 
point oniy, there remained the uncertainty of its distribution in the entire recess region, whereas, 
on the other hand, the joad could easily be set and held constant at a predetermined value. 


+ Flow rates were calculated fer a 0.055 inch diamet.. nozzle. The single nozzle was replaced hy twelve 


ung shea and thatr dismeters were cn 


: god fom S059 inches iv approximately U.uiow inches, 
tf As measured by a tap located at a radial distance of 6.344 inches from the center. 
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Yhe tests were conducted in the following manner, The load was applied to the bearing 
aah ine supply pressure adjusted tea prodetermined value. The recess depth was then increased 
ter h creased! quasiestaticaliy, until the beading becaiic unstable. In certain cases the procedure 
was reversed. eg. the recess depth was set and the sunpty pressure increased (or decreased) verv 
ateluadly, up te the threshold of instability. The axial oscillations were detected by tue primary 
eapteitince probe and viewed on the display of the cathode ray oscilloscope. Photographic rec- 
ards of amplitudes and froemencies of omcillatiuns were made for a number of runs. \ simuitane- 
aus check on previously determined. steady-state flow rates and gap widths was made for points 
it the immediate viciaity of the stability boundary. The bouadary was always approached by pro- 


ceeding from the stable region to the unstable region. 


‘The first. peeliminars stability test was performed, using a plane-faced pi -ton and a cen- 
trally located, 0.055 inch diameter supply nozzle. Results of this test are shown in Figure 4. 
The abrup: transition in the lower branch of the cucve, associated with the reduction in recess 
depth below the “cutoff” point (see Figure 37), provided yet another, valid reason for modifving 
the piston and, finally, replacing the single nozzle feed by a mu!tiple-nozzle supply system. ‘The 
upper Inaneh of the stability boundary constituted a marked departure from the predicted, single- 


Figure Hs contains results obtained with a modified piston, having a secondary, fixed vol- 
ume recess, but retaining a single, centrally-located supoly nozzle. It can be seen that the un- 
sttble region caclosed by the curve corresponding to a constant load test was smaller than the 
unsteble region bounded by the locus obtained when the load was adjusted in order to maintain a 
constant pressare at a tap, lucated 0.334 inches from the recess center. Both curves were marked 
by the absence of the sharp transition observed in the case of the plane-facer pistont, but contiay 
ued to display two characteristic branches, separated by the point at which their slopes became 
infinite. , 

\ secondary stability experiment was conducted in parallel with previously descrihed. 
static tests, the objective of the latter being to assess the reliability of the assumption that the 
recess pressure at constant load was not appreciably affected by changes in the supply pressure. 
It was observed that when the 0.0135 inch diameter, 0.125 inch long pressure tap (see Figure 
10) cnuneeted the recenm log relatively large, closed volumett, an effect, which could justly be 
described az hysteresis. was present. With the bearing stable at a particular setting of supply 
pressure aud recess depth, the supply pressure was increased very slowly and gradually, particu- 
lac care being exercised to avoid disturbances which might be large in comparisun with those ran- 
domly produced and always present in the fluid, foundation and the surrounding atmosphere. Under 
those conditions it was possible to increase the supply pressure considerably above the value at 
which the bearing became unstable if jarred, displaced or subjected toa sharp impact in the axial 
direction. The supply pressure was then reduced to a level at which the bearing became stable 
again, At this level, instability could no longer be induced, regardless of rate and magnitude of 
the disturbanve. Identical effects could be observed with the supply pressure held constant and 
the recess deptli being varied quasi-statically. With an isolated recess volume (e.g., closed pres- 
sure tap), no hysteresis effects were present, at least no measurable difference in limiting values 


Tt is to be noted that in Figure 45 the ordinate is the mean recess depth, that is the recess volume divided 
by the recess area, 


tt Consisting of drilled passages, boies and the volume of the hourdon tube. 
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ef stability pacureters was gute Result. of observations of this rather remarkable phenomenon 
hes Sits nee 
that the region Vanstable if disturbed” is quite 


tee shown in be tagure dan eis realized, af cons 


bhatt. the svstem could remain in the state of dastable equilibrium if all dis- 


rneoesy ne matter tea s calf, were absent. 


Tn the conduct of stebitity experiments described up to this point, the eniphasis was rather 
ono pLditative thea on qucttitative aspects of bearing performance. Subsequent tests, although 


: a. Sd ao _ 5 1 : oe Bet ‘ 
thea descrintion and te sults are presentedt in this section under the heading. of prelindaary ry 


were carried cut woth considerably ceeater accuracy and precision of measurements. Closer spac 

ine of intervals at which date was recorded and the determination of frequencies of anto-oscillas 

tions on the threshold of instability. contributed to the detection aud eluninaiion of certain troubie- 

some characteristics, inherent in the test tig, and subsequent improvements in correlation between : 


peedicted and experimental results. ; : 


Mability boundaries, obtained with a twelve supply nozzie, plane-faced piston, under con- 
stant load C472 1b) conditions, for 0.477 ~ 102 and 11,68 x 107? Ibesec? in effective bearing 
nidssest, are presented and compared with predicted results in Figure 1. Measured and predicted 
critical frequencies, ea. frequencies, 6. associated with limiting values of parameters on the 
observed threshold of instability an€ their calculated equivaleats, for which the roots of the char- 
acteristic equation QEb) on the axis of imaginaries were S - tie cr are presented and compared 
in Figure $8. Resalts afia test conducted under reduced load (222 lb) conditions are shown in 
Misuee $9. The stebility betndaries referred to in the foregoing, although similar in shape tu those 
shown in kigures th Ke and i, evhibited vet another, unusual property in regions denoted by 
‘resonance’ on all pertinent figures. “This peculiar trend in all loci, characterized by inflection 
points and abrupt changes in curvature, occurred invariably at frequencies between 220 and 225 
cps. indicating the existence of resonance in the system. 


The phenemenon described in the foregoing, a major source of discrepancy between theory 
and experiment, might wel! have gone unnoticed if critical values of supply pressure and recess 
depth hel bee recoraed at more widely spaced mtervals., Fortunately also, frequencies were re- 
corded in paraliel with other measurements, so that the local, abrupt transitions in several stabil- 
ity lect could be identified wich the s.une salue of approximately 220-225 eps. The preliminary 
dvaanic tests indicated a considerable acea ef qualitative and quantitative disagreement between 
predicted and eyperimental romalts, far in excess of what might have been anticipated on the basis 
of comparison between calculated and actual, static bearing characteristics, ‘There remained the 
problem of detecting possible sources of trouble inherent in the experimental apparatus and the ap- 
plicat’ a of corrective measures before final experiments could be conducted, Specifically, it was 
esoential te (a) identify: parts and components of the test rig assembly responsible for the exist- 
ence of resonance at particular frequencies of auto-oscillation and (b) to investigate possible in- 
teractioas between the motion of the bearing and tied flow in regulats, supply lines, 
ings, the loading estiader, etc. 


guide bear. 


ce 


a me cent tn t 


t The effective bearing roe consisted of the shaft and lower thrust plate assembly (Aporovimately 25 [b 
and $5 Ib weights respectively), f 
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3. Concluding Stobility Eaperinents 


This section begins wrth the description of preparatory cest.and adjustments of the experi- 
ment. apparatus which preceded and were of considerable impoctance in the final, experimental 
phase of the investigation. It contains a detailed report of concluding experiments and ends with 
an account of tests, performed as a check of repeatability of results. 


ft was suspected that the singular tends ass ciated with frequencies of auto-oscillations 
in the neighburhood of 220.295 cps, were caused by resununce of the entire test stand, considered 
as one large mass (Approximately 750 lb weight), supported on an elastic foundation, ‘The possi- 
bility of exciting one of the graver modes of vibrations in the massive, upper thrust plate was 
highly unlikely, bet could cot be entirely excluded, because of the rather indefinite manner of sup- 


port on three columns and levelling screws. 


The test bearing assembly was bolted to a very sturdy and rigid tripod (see Figure 18 for 
example), mounted directly on an asphalt tile covered, concrete floor. Crude deflection measure- 
ments, using a dial indicator and the experimenter’s weight for the applied load, yielded a natural 
frequency for the test riz mass and the clastic foundation of approximately 230 cps, a value ex- 
tremely close to the narrow tand in which the inflection points and rapid changes in curvature of 
stability loci were observed. In order to measure the natural frequency of this mase-spring system, 
a» well as the response of the upper thrust plate, it was decided to utilize the floated shaft agua 
seismic mass and the primary capacitance probe as the sensing element of a vibrometer, provident- 
lv, though inadvertently incorporated in the test rig design. A very minute, axial load was applied 
to the shaft thrust plate and balanced by the pressure in the loading cylinder. This furnished the 
desired ‘soft’? suspension for the floated mass and a correspondingly low, damped natural fre- 
quency of less than I cps. The base of the tripod was then pressed down and released, while the 
response was recorded by means of the oscilloscope camera. As anticipated, the recorded fre- 
quency was approximately 220-225 cps, 


The cause of the observed singularities thus firmly established, the entire test rig struc- 
ture was mounted on isolators. A record, obtained in a manner similur to that described in the 
foregoing, showed that the natural frequency of the isolated structure was 15 cps. A further, ef- 
fective reduction of isolator stiffness, such as could he ohtained with a Julius suspension, for 
example, or a sizeable increase in mass, was not feasible because of physical limitations.. The 
response of the upper thrust plate was also of interest. The plate was strack a light blow near the 


_center with a small hammer, and the resulting trace of the cathode ray oscilloscope was recorded, 
"Phe frequency observed in this case was approximately 730 cps, twice as large as the maximum 


frequency of auto-osciliations encoumered in ell experiments, but considerably lower than the mag- 
nitude estimated in the design phase of the apparatus, 


Photographic records of frequency measurements described in the preceding paragraphs are 
presented in Figure 50, in which ! representa the trace obtained with the structure resting on the 
floor, His the response <f the upper thrust plate and II refers to the structure when mounted on 
isolatws, ‘The uneven appearance of traces is due to the fact that, with relatively low pressures 
applied at the extremities, the shaft drifted cu: ‘inuvuoly in the axial direction, producing corre~ 
spunding changen in the electrostatic field hetween the prohe and the “seismic mass”. The os- 
cillations of the upper thrust plate, Il, are superimposed on the trace produced by vibrations of the 
test rig structure on isolators, both induced by the light hammer blow. 
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‘the concluding stability experiments were conducted with an isolated structure, so that »~ 
the zatlos uf lowest. observed frequencies of autu-oscillations and the natural frequency of the as- ; 


sembly-isolatur system were in the order of 14:1. Rearing alignment, lcakproofness, ‘zero settings 
of lustruments ete. were thersughty checked. \ tow frequency generator, the triangular-‘wave! 
output of which was used as u marker on photographic records. was checked against an electronic, 
digital counter and found to be stable and accurate within i percent in the 500-1200 cps band : 
width. Particularly great care was taken in recording readings of the recess depth dial indicator 
and of pressure gauges. Pressure readings were corrected to the nearest 0.05 psi for head loss 
between the measuring station and the point of admissien to the reservoir, upstream of the supply 
noceles. In approaching the critical recess depth on the threshold of instability, the piston was‘ 


lucked, so that final setting was made‘ by very gradually cvercoming a substantial amount of fric- - 
tional resistance. in addition to that offered by four sliding “O-ring seals, Whenever the stabil- oh 9 
ity boundary was approached by varying the supply pressure, the regulator control was. adjusted in | 


a manner, for which, perhaps, the most apt description would be “creeping”. 


It will be recalled that the stability analysis was haged on linearized equations of motion, 
implying among other things that. in the course of the experimental part of this investigation, 
seratbly small amplitudes of oscillation, possibly one percent of the equilibrium gap widtk,.would 
have ta be observed. Now, since sources of smull, random disturbances were always present in| 
the fluid, the foundation and the surrounding atmosphere, more or less well defined wave forme of: 
bearing response.-with amplitudes of the order of microinches, were always observed on the dis- - 
play of the cathode ray oscilloscope. On the other hand, no significant growth of oscillation was 
noted for very large chenges of supply pressure and recess depth, until their values approacied a” 
region in which relatively minute changes in either parameter were accompanied by an abrupt in- 
crease in the amplitude of vibration. In this narrow transition region the amplitude of vibration 


was modulated by irregular “beats”. Their incideuce, relative to the frequency of oscillation, oc-’. 


curred at more frequent intervals for large values of the critical recess depth. 


Perfect discontinuities and ideally defined boundaries seldom occur in physical problems. 
Nevertheless, the transition in ariplitudes from the microinch noise level to magnitudes. approach- 
ing the gap width itself was invariably accomplished by changes smaller than one percent of the 


. stability parameters designated as “‘critical’’. Thus the stability boundary is sharply defined, | 


Otherwise stated, the amplitude of oscillation remained very small and its rate of change (slope) 


with recess depth, ai constant supply pressuret, remained essentially zero, until both changed al= _ 
most discontinuously in a narrow region, the amplitude becoming very large and its rate of change © 


{elepe) almost infinite. 


The concepts and idcas contained in the preceding paragraphs can be summarized by stat- 
ing (a) that amplitudes of the order of 10°° inches were small compared with the bearing gap 
width, (b) tant the available instruments were capable of detecting and measuring amplitudes of 
the order of 10°> inches and, moreover, (c) that the order of magnitude of 10° inches was not - 
very critical as a critcrion of limiting values of recess depth and supply pressure at the inception 
of instability. oe 


oe 
TAltematively, with supply Pressure at constant recess depth. 
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The fouregeing statement may he substantiated by referring to Figures 51, 52, 53 and 34, 
which i!lustrate the abrup: transition in amplitude growth when approaching the boundary from the 
atable side. by increasing (ur decreasing) the recess denth at several settings of tt.: supply pres- 
sure. The lower photographs in each row were obtained with long exposure times and slow sweep 
rues, They show quite clearls the irregular “beats” referred to previously. The upper photo- 
eraphs display wave forms of auto-oscillations in random time intervals and were obtained with 
sweep cates one thousand times faster than corresponding photographs in the lower rows. 


Stability boundaries and frequency determinations, obtained in the course of concluding ex- 


: wriments, are presented and compared with predicted results in Figures 35 and 56. A correspond- 
H i ft P 


ing tabulation of experimental data is included in Appendix 4, following the table of calculated 
results, Photographic records of frequency measucments are shown in Figures ST and 58. In the 
presence of irregular beats, described in the foregoing paragraphs, particularly at large values of 
the recess depih, the determination of critical parameters was unavoidably subject to judgment of 
the observer, The procedure fullowed by the experimenter was to obyerve the truce on the cathode 
ray oscilloscope while adjusting the recess depth, without noting the reading of the recess. depth 
dial indicator, The objective was to avoid bias in selecting the point ot which auto-oscillations 
were estimated to have reached amplitudes between approximately 1 x 10°5 to 2 x 10°5 inches. 

An alternative method was used for a number of points in which the stability boundary was ap- 
proached by varying the supply pressure instead of the recess depth, thus providing a check on 
the reliability of measurement. In the latter case, the attention of the experimenter was focusned 
on the oscilloscone display without observing the pressure gauge reading until the critical point 
was determined. No significaat differences were noted, regardless of the manner in which the sta- 
bility boundary was approached, indicating that results were reproducible according to the select- 
ed criterion, as well as rather insensitive to inexactitudes involved in its selection. Several 
measurements were repeated un days when the barometric pressure and the ambient (room) temper- 
ature were slightly higher, or lower, but the observed differences were negligibly small. It is -~- 
timated that the repeatability of limiting stability parameter measurements is at least within L to 2 
percent of readings obtained in the course of conelnding experimental runs. 


The very substantial and, in this cane, rather encouraging qualitative and quantitative dif- 
ferences in results of preliminary and concluding experiments were attributed éatirely tu the igola- 
tion of the test rig structure. since no other changes were effected in the system. On the other 
hand, all static and dy namic experiments seemed to indicaie that, under certain conditions, even 
minute changes in system characteristics had an enormous effect on the stability of the bearing. 
In the final phase of this investigation, therefore, the experiments were accomnanied by. numerous 
and important secondary tests, with the objective of (a) discovering posible causes of remaining 
divergencies hetween theory and experiment, particulas!; xuch as inigh. have accounted for the 
double-valuedness of stability loci, or, failing this, (b) adding weight to the reliability of the con- 
cluding experiments, Those secondary tests will now be described, 


Ail experimentally determined stability boundaries were double-valued and exhibited an 
“‘upper’’ branch, not predicted by the present analysis. Although no ca'culated data were avail- 
able for recess pressures other than 73.5 psia (corresponding to a 472 lb load), a stability test 
under reduced load conditions (222 Ib), cursar 'y mentioned in a preceding section, wan carried 
out. The reduned laud toot war dictated by the demisability of ascertaining whether the double. 
valuedness would also persist in the case of a locus corresponding to lewer Reynolds number, 
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ures Wand 49, in conjunction with reference te Figures 41 and 42, clearly indicaies the presence 
of “upper branches’ in loci corresponding te different ranges of.the Reynolds number. 


Next. effects of misalignment were investigated. This investigation was accomplished by 
two methods, First, the constraint of the shaftt was removed and the relative positions of the up- 
perand lower Ususi plaice were anered, These changes prodiced a very substantial -nisalign- 
ment duc to the run-out of the lower thrust surface with respect to the shaft joumals, Although 
ver minor differences in limiting values of supply pressure and recess depth were noted, no sig- 
nificant change in the trend of the stability boundary was observed. Similar conditions prevailed 
when the rotor was brought up to a speed of approximately 2500 rpm, by means of an air jet direct- 
ed onto the turbine wheel buckets. In the second misalignment test, the supply pressures of the 
upper and lower guide bearings was varied in such a wanner as to tilt the shaft within the limits 
of tearing cloamnce. In addition to wisalignmer.t, the accompanying effect, in this case, was to 
change the ‘stiffness’? of the guide bearings. \o substantial changes in stability characteristics 
were observed, except when the supply pressure to one, or both bearings, was reduced to the ex- 
tent which precluded their functioning as effective guides and constraints, or when metal to metal 
contact between the bearings and journals was established. 


Experiments conducted with the isolated and non-isolated test rig assembly not only dem- 
onstrated local effects in the vicinity of the resonance point, but also the presence of large differ- 
ences in limiting values of stability narametere for freancucies uf incipient aute-uscillation von. 
siderably removed from the critical. The existence of other critical frequencies in the experimen- 
tal band width and the effect cf smal! vibrations of the upper thrust plate could not be entirely dis- 
missed. It was, therefore, decided to stiffen the upper thrust plate. This was accomplished by 
placing three 3 4 inch diameter jacking screws between the column-bracing, upper plate and the 
Upper thrust plate (see assembly drawing Figure 7) in Positions straddling the columns and as 
close to the bearing center as was feasible. The jacking resulted in a slight amount of controlled, 
symmetrical deflection. It was known, however, at-that time that no significant changes in stahil- 
ity characteristics were associated with a moderate degree of misalignment. The rigidity of the 
upper thrust plate thus considerably increased, measurements for a series of previously estub- 
lished, experimental points were repeated, without, however, any appreciable departures having 
been noted, 


In the analysis, as well as in the course of experimental investigat. 4s, it was postulated 
that the load, and consequently the pressure in the loading cylinder, remained co:ustant at all 
times. Accordingly, the. design of the loading cell was such that fractional press ire changes, as- 
sociated with motion of the shaft, were to be of the order of smallness of the displaced to total 
evlindcr volume ratio. (Approximately 10°°.) There existed the possibility, how: er, that the 
lower end of the shaft, acting as a piston type source, gave rise to a complex, ‘stationary wave 
system. If wave reinforcement occurred, the effect could be appreciatle and divergence between 


the theoretical and the physical models signifi-ant. The possible presence and influence of res- 


onance in the chamber couid he confirmed or disproved hy aleveing the shape aud volume of the 


Tt Two jets impinging on each side a guide fin, mounted on the vhaft collar. See Fignee 18 
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Na 


lading ov tinder, Vsieple method was employed. in owhich the cy Linder was fitted with as 
siaudpipe and ther filled vith water in increments up to three quarters of the total volume, 
signitionue departure se however, fren previously established results were observed. ) 


The valitits of ine concept ofa “constant pressure supply reserveie” depends on the ab- 
sence of pulsations transnitted through the supply restrictors, The likelihand of propagation aud, 
possibly. of c-fiection of pressure disturbances apsirean of the supply nozzles could not be dis- 
missed. Unfortunateiv. insielhetin: -f saitably siall and sensitive pressure pick-ups, or hot wire 
anemometers, in the immediate vicinity of the nuzzte eptrance and exit would have involved sub- 
stantial changes in design ad elaborate. additional insteementation, Subsequent investigations, 
therefore. were Limited te a series of simple checks, intended ta locate possible sources of trouble 
inthe air supply avstem. [tis. perhaps, both useful and instructive to render a brief account of 


sone of the steps taken in this “trot ble-shuoting™ procedure: 


a. Vitwo cubic foot plenum chamber was installed at a point where the supply line was 
joined tu the upper threst plate (see assembly d:awing in Figure 7), in order to increase 


the capacity of the supply reser: nd provide a zone ip which small fluid disturb- 


vaces couid be dissipated. 


he The compressor was shut off se that no unsteadiness uf flow originated at the source 
of the supply while readings were recorded. The capacity of the main supply reservoir 
was sufficiently large to maintain a constant supply pressure in a time interval requiced 


te perform one of several measurements a. 


© Since instabilities of pressure regulators and, ppeumutic valves are similar ‘in nature to 
hose occurring in externally pressurized bearings, the regulators and valves, placed in 
series in lines supplying the bearings and leading cylinder (see Figure 17), were ad- 
justed for varying degrees of pressure drop across the metering elemeute. Various set- 
tings were made and a number of measurements were repeated. 


None of the foregoing measures had any effect on the stability characteristics of the test 
hearing, : ‘ 

Valid comparisons between theoretical and experimental results could only be made if the 
pitching (cocking} motion uf the shaft vemained smaii in comparison with the amplitudes of axial 
oscillation. In the course of several preliminary and final runs, the auxiliary capacitance probe 
was mounted in different positions arcuad and along the shaft. but no measucable amplitudes were I 
detected. Vreliminary calculations curried out in the design phase of the apparatus, in which the 
shaft and the guide bearings were consider-d as 4 simple, two-degree-of-freedam system, indicated ‘ 

‘that the slower mole of vibration could possibly be excited by frequencies of auto-oscillation in | 
the upper range of the experimental band width. Figuie 39 shows results of a vibration test, in 
‘which the auxiliary capacitance probe was located in a positiont where large an.plitudes at tre- | 
quencies close to that of the fiist mode could mont readily be observed, Since the sensitivity of ' 
the auxiliary probe was considerably less than that of the primary probe, the gain of the primary | 
instrument was decreased and the amplitude of auto-oscillations increased to 4x 1074 inches, ap- 
provimately twenty times the value observed on the threshold of instability. This permitted a ber- 
ler comparison to he drawr between amplitudes or axin| ond pitching motions. The frequency af i ‘ 


t Along the six inch outer diameter of the shaft thrust plate, 
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aulowosciliations was then varied fram 375 to 200 cps, The second photograph in Figure 59 shows 


the largest amplitude of pitching which was observed, occurring at 2091 rad/sec (333 cps), It is’ 
that more than two modes could be excited, that the shaft could perform small, 


tealtped of CPE 

conical gveations within the limits of bearing clearance, or momentarily rotate and translate, put 
net necessarily im the piane of the capacitance probe. Nevertheless, Figure 59 shows the largest 
froximately 100:1) ever observed in any of the monitoring positions and there is 


amplitude ratio (ap 
at the overall effect of pitching was a very minor one. 


little doubt. therefore. th 
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Vi. DISCUSSION 


1. Comparison of \nalvtioul and by perimental Results 
! 


Comparison of thearetical anl experimental loci in Figure 55 indicates ronwoaniely good 
detcement in the regen in wich pressure differences between the supply reservoir and the recess 
are highest. In this recian, values.of 8 obtained from the present analysis represent a consid-| , 
ereble improvement when compared: with those predicted from the simplified, lumped-parameter ' 
theory, (See Ligure 4.) Ys the pressure drop across the supply nezzle decreases, the deviation 
of empirical curves from their cafenbued equivalents becomes more pronounced. Beyond the 
points of vertical tangency neither quantitative nor qualitative agreement exists. 


With regard to slope, empirical curves and those predicted from the present analysis agree 
fairly well over a considerable part of te-experimental range. The same situation does not hold 
true forthe single lacus obtained from the simplified theory. The present analysis: predicts 
conservative (safe) limits, whereas the simplified analysis gives values which are either too high 
ot tos low. Atmost significant improvement, introduced through the refinements incorporated in 
the distributed-parameter approach, is that the present analysis gives a correct interpretation with 
regard to influence of the effective hearing mass on stability, a fact to which experimental results 
lend unequivocal support. As emphasized previously, one of the major shortcomings of lumped- 
parameter analyses is that they predict the same locus, regardless of magnitude of the beacing 


t 
Mass, 


It is interesting to note that forthe larger mass, agreement between predicted and experi- 
mental loci in Figure 55 is better and extends over a larger portion of the curves. It is more than 
nrobable'that, as the bearing mass increases, agreement between predicted and actual results pro- 
uressively improves. This inference may be particularly important, because in actual applications 
one might expect considerably higher than 0.65 ui 0,10 weight to thrust ratios. Since from a prac: 
tical point of vew one is usually interested in “safe depths” of the order of 5 x 107? to 20 x 1074 
inches, it is encouraging to note that is is precisely in this region where'predictions based on the 
present analysis constitute a considerable improvement in comparison with those based on simpli- 


fied analyses. 


2. Remarks on the Relative Validity sf Theoretical Assumptions and Their Applicability 
To the Experimental Model 


The present investigation included both, analysis and experiment. To the investigators, 
who have engaged in this daal undertaking, falls the difficult task of accounting for disagreement 
in areas where it exists. In what follows,.an attempt is made to inquire into likely causes of 
divergence between analytical predictions and empirical results, particularly the single and 
double valuedness of the respective stability loci. 


It will be recalled that the dynamic lubrication equation (4) was required to satisfy a 
boundary condition at the outer bearing perimete: .nd to be compatible with the continuity equation 


T The bearing mass becomes @ significant stability parameter in presence of appreciable, extemal damping 
only. (See Appendix 2.) 
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(lof the variable-valume recess. Now, in implementing the continuity equation, it was assunied 
that he infl.s throneh the supply nozzle. W. and the outflow along the perimeter of the recess, 
Sy. coubd he considercd as quasi-steady, The quasi-steady approach has been successfully em- 
plowed in proldens of charging and discharging of fixed and variabie-volume eslinders, in research 
related te the dynamics of pneumatic controls. A typical example is the work of Shearer [9!. The 
flow throuzh orifices, under pulsating downstream conditions, has been investigated, both theoret- 
ical. aad esperimenially, by Schultz-Gremew [16 He sugeested that limits of applicability of 
steady -Hoaw equations could be assessed from the reciprocal of the Strouhal number: 


Ee 
g-— 

qd 
T = Average velocity 
t, - Period of pulsation 
d= Orifice diameter 


which is chrectls related in unidimensional, inviscid flow to the ratio of local and convective 


acceleration components, Us vn 'fe concluded, that jor values of & > 510, the steady- 
ie f 


flow equetions vielded accurate values of instantaneous flow rates aed the effect of local arcel- 
crution was negligible. \ithough maximum values of & obtained in the course of present experi- 
ments were larger that 10%, magnitudes of @ were generally comparable with 510, and thus not 
sufficiently far removed froma transition zone, in which the effect of pulsation may have been 
siunificdnt. Other investigators discovered that pulsations propagated upstream of the restrictor 


and gave rise to complex wave patterns and phase relations on either side of the orifice [47, 48, 
19,59], ; 


The adequacy of considering infiow as quasi-steady is questionable, at least in certain 
tanges of pressure, frequency and geometrical parameters of the system. In,the limiting case, when 
PP approaches the critical ratio, no pressure pulses propagate apstream of the nozzle throat. 
li ane assumes, however. that acoustic equations hold when P_/P. approaches unity, the theory 
of horns, to wich the snpoly nozzles bear a geometrical resemblance, predicts definite phase 
relutions between pressure and flow at the throat [64]. The existence, therefore, of appreciable 
time lags cannot be entire!y dismissed, in which case the as~:med, “‘in phase’ perturbation : 
relation, w, = Cgp,, may be deficient, 


If we now direct our attention to the outflow from the recess, the entrance efiects, associated 
with fluid acceleration and flow separation at the sharp edge of the perimeter, sie further com- 
plicated by the variation of the area of admission to the slot under dynamic conditions. [Here, 


the use of equation (9), e.g., W, = - (constant) Wns. implied not only that the outflow could 
r 


he considered as quasi-steady, but also that there was no transition and that a viscous, laminar 
flow regime became at once established as ti.e fluid entered the slot. Consider also the case of 


small-amplitade vibration when the recess volume is large and the pressure deviation, Pee tenda 


to become very small in comparison with p, the perturbation pressure in the slot. In fact, due to 
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the trapping and squeezing of the fluid in the slot, maxima of p are likely to occur close to the 
fecens perimeter, where the gas is relatively dense. ‘This may cause large variations of the 
perturhation pressure gradient in the immediate vicinity of the recess and, in turn, induce appre- 
ciable fluid acceleration and deceleration in the region of admission to the slot.t Moreover, the 


3 - 
o 


Hnearized expression, Wa = Ch - Cop - Co. corresponding to equation (9), does not take 
Be i ‘ z 


into account the complex phenomena occurring in the neighborhood of the recess boundary. 
Divergence of results may be due partly te che Cumulative discrepancy involved in expressing the 
difference, w) - wy, corresponding to the nghthand side of the continuity equation (7). 


Next, we consider possible shortcomings of assumptions involved in the mass accumulation 
teem onthe left-hand side of the continuity equation (7). If we now focus our attention on the inter- 
nal peveesses occurring within the variable-volume recess, such as the impingement, deflection or. 
disintegration of jets and the interaction between decelerating and accelerating fluid streams, 
the physical picture defies mathematical description. Unavoidably, therefore, it becomes neces- 
sary to resort to simplifying assumptions, or postulates, sack as unifonnity of ‘“average’’ prop- 
erties at any instant of time. The latter assumption was made in the qualitative analysis of 
periodic charging and discharging of the variable-volume recess in Appendix 5. 


Granted the high decree of oversimplification in formulating the problem, the analysis of Ap- 
pendix Snevertheless suggests that the pressure and density in the recess ure not in phase. Numer- 
ical examples, indicate, however, that no appreciable time lag is involved, so that the process 
may he reasonably well described by P p™= constant. Although the exponent, n.is then not a 
function of time, it does depend on a dimensionless “charge” parameter, Q = W h.V;o% in 
which Wo. p+ ¥,, refer to the equilibrium flow rate, the recess gas density and the recess 
volume, and ¢ is the circular frequency of oscillation, 


Considerable numerical discrepancies between theoretical and experimental results could be 
accounted for by selecting mare “representative’’ values of the exponent n, but it is by no means 
implied here that an explanation for the double valuedness of experimental loci could be supplied 
by assigning to n.values ather than y = I. tin the recess, or ather than n ~ Lin the slot.tt 


It is considered that the inertia forces in the slot were small in comparison with viscous 
forces. Although no definite, analytical proof has been furnished to support this contention, the 
derivation in Appendix 6 suggests that the neglect of inertia contributions in the hearing clearance 
may not be verv serious, and that the dynamic lubrication equation (4)is a good first approximation, 
Moreover, in the presenve of significant inertia effects, one would expect the disagreement be- 
tween predicted and empirical results to be greatest at relutively high frequencies of oscillation, 
whereus exactly the opposite is indicated by the experimental evidence. 


{In the foregoing discussion, an attempt was made to pmbe into the more questionable aspects 
of the present analvsis. A critical inquiry would he incomplete and one-sided, if one omitted ta 


t An analogo1s condition is illustrated diagrammatic Wly in Figure 6a., showing a gre 
instantanecus pressure profile and positive gradient at the edge of the recess, 


tt in the slot, P/p = constant was the acsumed relation. See also Appendix 6, 
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examine both the experimental-model and the procedure for possible deficiencies und inherent 

‘ : — ‘Con. 
sources af error, Here, however, che numerous, »«condary tests, described in the section Con 
clading Stabjlitv Faperiments.” have narrowed down considerably the are2 af search for a feas- 


ible explanation, 


The turning point in the experimental phase of this investigation was the discovery of reso- 


nance in the avstem, [sofation of the test structure nut only altered the shape of stability loci near 


altered the 
fesonauce points, but alse displaced the curves in a direction effecting considerably better agrce- 
mer 


between theory and exper'ment. The extent to which actual motion departs from the ideali- 
m wle-degree-of-freedom system, for which perfect isolation, absolute 
Ligudity aad absence of external damping are assumed, is most difficult to assess. One need 
z that 10° © inches represented 10 percent of what could be meaningfully inter- 


preted in the presence of mechanical and electronic noise. Since the capacitance gauge meas- 


Zation of 


undamined, st 


ured the relative, instantane: 


displacement of bearing surfaces at a point, it provided no 
information with regard to the absolute motion of boundaries. Nor was it feasible tu furnish a 
perfectly isolated and immobile reference frame for external mounting of capacitance probes. 


Thr extension of the simplified analysis of Appendix 2 to a damped, two-deprees-of-freedom 
avstem suggested the possibility of a double valued stability locus. An approximate value of the 
critical damping ratio was estimated from Figure 50.f Calculations then indicated limiting values 
of critical recess depth on the upper branch of the locus to he of the order 107 inches, as com- 
pared with 107! inches in the case of experimental results. The simplified analysis. of course, 
lends oaly qualitative support to the existence of double valued loci and may predict grossly 
exagecrated results in this region. It is not unlikely that the extension of the present unalvsis, 
to include the effects of elastic and damping properties of the entire system, would result in 
moce realistic predictions of critical values of stability parameters. \t any rate, one cannot rule 
out that the divergence between theory and experiment was due to appreciable motion and accom- 
panving dissipative effects in the supporting structure, whereas perfect rigidity and notion of 
one boundary only were hypothesized for the idealized modet. The importance of considering the 
motion of the supporting structure was stressed ir ~ recent naner by Tsai [52], who investigated 
the stability of a simple, pneumatic pressure reducer. The author states in the introduction: 
Also, the vibrations of the mounting structure may induce undesirable oscillations in the re- 
ducer. To obtain a complete picture of performance, one must therefore investigate the dynamic 
characteristics of the reducer and analyze the problems of its natural frequency and stability.” 


Of equal importance is the extent to which the actual, steady-flow pattern conforms to that 
which has been postulated for the idealized model. The flow of compressible fluids in narrow chan- 
nels and in externally-pressurized bearings, in varivuus Reynolds and Mach number regimes, has 
heenthe subject of numerous, theoretical and experimental investigations.tt The Re-nolds and ach 
numbers corresponding to the flow range of the present investigation were sensibls stmsll, and 
deviations from calculated, steady-flow bearing characteristics were of the order 2-5 percent. 
Nevertheless, once oscillation is started, the influence of entrance phenomena on stability may be- 
come of far greater significance than their effect under static conditions would seem to indicate. 


T Trace of decaying vibration of structure when mo ~ted on isolators. 


TT Representative examples may he found in references [19, 20, 21, 22, 23, 24, 27, 31]. A rather complete — 
list is contained in a bibliography by E.R. Sciulli, reference [44]. 
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The flow pattern in the entrance region does sei couform to the idealization of the experi- 
mental medet. te the flaw ground the sharp eae of the recess, separation takes place and the Tnid 
rea converses, a toroidal vertex forming beaeath the contraction, \cGinn (35! ohtaiaed stecam- 
line patterns for water flaw ina sharpeedued, radial diffuser by means of dve filaments, and re- 
corded pholugraphroadits the vortices formed downstream of the separation point. Shires [225 
using air, demonstrated the existence of pressure depressions in the entrance zones of several. 
pacailel-surface slot configurations, for various pressure ratios and slat width. The Revnolds 
number in the experiments of WeGinn and Shires was af the same order as that in the present 


investigaulion, 


Consider a vibration, induced bs a random disturbance and accompanied by changes of pres- 
sure gradient ip the entrance region. The vibration will cause alternate elongation and contractions 
of the transverse dimension of the vortes ring. The outflow from the recess depends an the width 
of the effective fiew area, rather than on the width of the slot. The amplitude and phase relation 
Letwoen the slot width and the width of the contraction may not be a simple one, and the repre- 
sp 
rf muy lead to 


sentation of the instantaacous flow rate by the expression Wy = —(canstant) 3p. 
f } r 

considerable error Conceivably alse. an appreciable amount of energy, which would otherwise be 

dissipated in the process of a self-sustained oscillation, may instead he absorbed through en- 


irunee losses. The entrance zone mav thus provide an internal damping mechanism and have a 


stabilizing effect on the hearing. Siace entrance phenomena were net considered in the analysis, 


theory may predict instability in regions which experiments show to be stable. 


On the other hand. as one follows the lower branch of a lacus and departure fram predicted 
results increases, the Revnolds and Mach aunbercs d. crease, so that entrance effects became less 
pronounced, up to the point of vertical tanzenes of the locus. Thereafter, the upposite holds true, 
It would be inconsistent, therefore, to attribute the main cause of divergence of results solely to 
entrance effects, except to sav that their discegard may have heen a contributing factor. 


In this section. a number of likely, but by ca means proven sources of divergence between 


theoretical and empirical results, hve been discussed, \though the observed discrepancies may 


have been a cumulative result of seveai deficiencies involved in the theoretical assumptions, a» 


well as those inherent in the experimental apparatus. the most probable sources of errar which 
suggest themselves are: 


a The nawethle te cafficience ifthe quasi-steady flow representatioa and of the pustulation 
of uniform propertics in expressing the continuity equation of the variuble-valume re- 


cess, 


b, Aa appresiable effect of the relative motion of the bearing surfaces, associated with 
imperfect isolation and or the distributed, clastic and damping pmperties of the sup- 
porting structure, 


c. The presence of entrance effects in the experiment«] model, contrasted vei neglect 
in the present analysis, : 


AW 
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VEE. CONCLUDING REMARKS AND “RECOWMENDATION NS 


It has been demonstrated that the dynamic lubrication equation (4) gives a more faithful 
representation of the effects of partial trapping and squeezing of a compressible fluid film in a 
long and narrow stot than could hitherto be obtained with the aid of arbitrary assumptions with 
regard ta the flaw and the pressure profile. Moreover, the present results have been obtained 
without resnet ta any corre tion factors, It has alse, been shown, both theoretically and experi- 
mentally, that the bearing mass has a significant effect on stability and that stability is en- 
hanced by decreasing the bearing mans. As a consequence of results of the analvtical and the 
experimental phases of this investigation, it has been possible to assess more critically the 


earlier, lumped-parameter theories, 


Experience gained in the course of experiments has shown that the response of the bearing 
can he greatly influenced by the method of mounting, indicating that in actual applications the 
mass aad the elastic and damping properties of the supporting structure may have to be considered 
in tormulating the stability problem. Stability has been shown to be extremely sensitive to minute 
changes uf certain hearing parameters und svstem characteristics and quite insensitive to others, 
‘The experimental results have also substantiated the validitv of the linearized approach in the 
analysis, since, in regions designated as “‘stable,’’ no self-sustaining oscillation could be in- 
duced, recardleas of the magnitude of disturbance. 


There are several areas in which the present research could be further extended and sup- 
plemented. The stability of bearings cannot be treated as a separate problem and must be dealt 
with in conjunction with other aspects of an entire svstem. [t is not always possible, however, 
to formulate the equations of motion of a complex system in terms of distributed parameters. 

Even after such formulation, it may be impossible to obtain solutions. Moreover, system consider- 
ations aside, the distributed parameter approach may also become too involved for more complex 
bearing configurations. Clearly, therefore, a mare rational method uf ieohea paremotet represen- 
tation would be desirable. 


The present analysis suggests an approach in this direction which merits further consider- 
ation. It will be recalled in the case of the rectangular geometry, that in the series representa- 
tion of the functions Yy and yy] we capitalized on the smallness of the parameter € and neglected 
higher order terms of the expansion. When the first four terms of the expansion are retained, the 
characteristic equation reduces to a polynomial of the third degree in s. On the other hand, it has 
been shown in Appendix 2 that the characteristic equation obtained from the simple analy sis is 
also a polynomial of the third degree in s. Since both equations are of the same form and differ 
only in their coefficients, the results of the present analysis may be helpful in devising a 
rational method of correlation of the distributed and the lumped-parameter representations. In the 
case of the circular bearing, it may also be feasinle to obtain solutions y, and yj, in the form of 
series expansions in terms of the small parameter, 12uR ayP. oll as which characterizes the 
equation (19b). [t may then he possible to represent the ‘characteristic equation (31b) by a paly- 
nomia! ins, following a procedure analogous to that which was employed in the case of the rec. 
tangular bearing configuration. A correlation cf coefficients in the characteristic polynomials of 
the present and rimnlified analyses wanld vive the advantare of farmetating the nehlem in terme 


of ordinary, rather than partial differential equations. 
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ft would be iegedie te retain the quasi-steady flow representation in furure stahilitv | ' 


Hieits of such an approach. Nevertheless, in certain flow 
cal description of the charging — dis- 


avalyses, because of the mbkerent sian 
WW become necessary to modify the mathems 
ess of sariable vofumes in bearings. Perhaps, as a first step in this direction, it 


ces of simple noseles aud of orifices, in a inanner 


ve 


“out 


charging proc 
woithd be useful to determine the throat itupe: 
analogous to derivations in the acoustic theory. A better understanding of the phenomena 
eceurring in the entrance region of a siot is also needed, before an attempt is made at an improved 
aithenatical description of fluid How in the enteance zone. A studv of streamline patterns and 

af vartes formations should be made, when the boundary oscillates in the transverse direction of 


the Slot. A complete time record of flaw pitieas conid be obtained with the aid of a high speed 


camera, using otherwise methods similar to those emploved vy McGinn [35] for the case of steady 


flow of water in the entrance zone of a sharp-edged, radial slot. 


Advance consideration has been given to future experiments which would involve rotation 
of the shaft. For this purpose, an air-driven turbine has been incorporated in the design of the 
shaft, details of shich are shown in Vigure 5. The ettect on stability and other rotor character- 
istics due to pressure cradients, induced through centrifugal forces and surface misalignment, is 

= Ae of cousiderable interest and should be studied. Also, additional stability loci could be obtained 
for various pressure ratios and types of supply restrictor, under rotating, or non-rotating con- 
‘ditions, “The experimental apparatus is capable of providing a great deal of valuable infarmation 
and can also be modified for the purpose of experimentation with other bearing configurations. {2 
ix recommended, therefore, that full use be made in the future of the cuisting facility. 


Minally, since the simple geometries considered in this investigation are frequently en- 
countered in practical appiications, and since god agreement exists in regions which are of in- 
terest, the programming of a digital computer, to provide design data for a representative range of 
parameters, is warranted and recommended. 
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SYMBOI 


{ ' 


A. 
A 
A, 
a 


(sa) 


1X. NOMENCLATURE 


MEANING 


_ “Cgnivaleut’’ arca, defined in Appendix 2 


Recess area 
Coefficient, - (P.2-P2)/X-X, 


Complex function, defined in text 
R 
Coefficient, - (P,2 - P,7)/log. Rr 
r 


Coefficients, defined in text 


Damping constant associated with mase M, 
Appendix 2 
Specific heat at constant pressure 
Specific heat at constant volume 
Nozzle throat diameter 
2.718... 


Base of antural logarithms, 2.73 


Forcing function 


Perturbation forcing function in text, (H(t)/H* in 


Appendix 6) 
Functions of a and P,, defined in text 
Bearing gap width 
Bearing can width, eqnilibriem 
Bearing gap width perturbation, H ~ H, 


Bearing gap width perturbation amplitude 


Pressure integral, defined in text and expressed 


in Appendix | 
v-1 
Bessel function 


Spring constant associated with mass M, 
Appendix 2 


Coefficients, defined in text 


Thermal conductivity 


FIRST 
APPEARANCE 


A2-1 
A2-1, AS—6 
16a 


14, (A6—7) 


_ 43b 


aod 
~ 


ww UNM 


1 
AS-6 
26b, Al-3 


15a 
AR_A 


a ee, 


_-1-A2049-12 


UNITS 


in? 


th?/in 


1b2/in* 


ib-aec/in 


in?/sec2=°F 
in?/sec2—°F 
in 
No dim. 
ib 
Ib, (No din.) 


wer we 
on we 8 
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sy inor MEANING saan UNITS 
H Bearing load in text and in Appendix 1(R - RB, Al-9,(A6—7) Ib, (in) 
> : in Appendis 6) 
< u thiseaC eaneacainpon A2-10 Ib-nen?/in | 
eee MI Total gas mass-content of bearing A2~4 lb-wec? ‘in i * 
M Gas maes- content of recess volume ; ; ? Ib-sec?/in 
m ; tearing mass (Mevable, or effective bearing 6 Ibesec?/in 
mass} 
4 Pressure in hina gan, abs, 1 Ib/in? 
of a Piss Atmospheric (ambient) pressure, ahr. . § lb/in? 
Pp _ Pressure in bearing gap, abs., equilibrium 6 Ib/in? 
Pp. Pressure in bearing recess, abs. 6. Ib/in? 
Pio Pressure in bearing recess, abs., equilibrium 6 : ib/in2 
1. Supply pressure, abs. 8 ~— Yb/in? 
. p Perturbation’ pressure in bearing gapt, P ~ P, 12 Ib/in? 
B, Perturbation pressure in bearing recess, P,/-P_ 12 lb/in? 
Q) Complex function, defined in test 2B ‘No dim. 
G “Charge’’ parameter, defined in Appendix 5 AS—11 No dim. wi 
R Radiuy of bearing Bb in 
KR. Radius of recess” 6b in 
a Gas constant, Cp7 ey 8 ; in?/sec7=°F 
r Radial ccordinate ih : in 
s Complex variable 19. L/sec 
T . Temperature, abs. A6--4 : oF : 
T, ‘Temperature in bearing gap. abs., equilibrium 8 oF | 
(T, = T,, <1, =T,) 
Tine 3 sec | 
tr Velocity functions, defined in Appendix 6 A6~15 No dim. 
a Component of velocity in the x and/or t-direction . A6~1 in/sec | 
. u Mean velocity in the x and/or redirection ; 1 ; in ‘sec | 
7 Eacept im appenax >, 1m which p = P,~.P,, is the perturbation pressure in the recess und the sabscript 
‘ **s' is omitted. 1 
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SYMEOI . IEANING , er tee UNITS 
- Velocity functions, defined in Appendiv 6 A615 No dim, 
‘a Reems value AS—6 in’ 
\.. Reccas welunaaiuiliiniem ASA ind 
Pressure variable, defined in text 18 
i Bearing gus mass flow rate, equilibrium 16 Ib-sec/in 
‘, Rate of gas inflow through supply nozzte , i "  |besec.’in 
Me Rate of gas outflew from recess for outflow from 7, A5—1,(A2-4) ——Ibesec/‘in 
bearing) ’ ' ; 
TW, Velocity functions defined i Apnendix 6 - AGH=15 ‘No. dim. 
ee Velocity component in the z-direction A6~1 in ‘sec 
w Perturbation inflow, v - K A5-6 thesec, in 
Way Pe turbation ouflow, 8, - wo ASG ib=sec/in 
\ Half width of rectangular bearing . 5a in 
\, Half width of rectangular recess 6a in 
\ Coordinate Mad la , in 
‘int Solutions of homogeneous part cf equation 19 28, 32a, 37b 
Substitution variable, z = P? 18 Ib? ‘in 
z Coordinate A6~1 in 
a Compressibility parameter, 1QyaR? /P_H,? 3th No dim. 
B Parameter, 1] - (7 VE : 37b, Al-3 No dim. 
P © Rr | 
B Parameter, H,/é + H, AS-6 No dim. 
y Ratio of specific heats. c,/c, a 8 No dim. | 
5 Recess depth 10 in 
&. Critical recess depth 44b ia | 
2 4 : 
€ Coefficrent, biel fx-x,\ , text a in®-sec/Ib? . : A 
n? \ p2-p2 J 
Ae Cee | 
€ Dimensionless amplitude, h,/H,. Appendix 5 A5-6 No dim. : 
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SYMBOL, _' MEANING eeerka AGE UNITS 
. Parameter, H°/K - R_. Apperdix 6 A6-7 No dim. 
¢ _ Dummy variable, corresponding to z = Pp. 28a Ib/in4 
n Coefficient, 12,/H2 : 2b Ib-sec/in 4 
A, Temperature functions, defined in Appendix 6 A6-15 No dim. 
Atos Coefficients, defined in Appendix 2, A2-4 
og Coefficients. defined in Appendix 2 A217 
i! Dynamic viscunity a i lb-sec/in? 
7 v Kinematic viscosity A6—7 in?/sec 
é Dimensionless radius, r/R, ; 37h No dim. 
paso ; n Constant, 3.141... 8 No dim. 
a Pressure functions, defined in Appendix 6 A6-15 No dim. 
Pp * Dummy variable, aicmandne tor . 2b in 
& a p Gas density in bearing gap a8 Ib-sec?/ia* 
oui P, Gas density in bearing recess 10 Ib-sec 2/in? 
Pro Gas density in bearing recess, equilibrium ll Ib-sec?/in4 
o Perturbation density, p, - Pro A5-6 Ib-sec?/in4 
g Suit Real component of Yun 38b 
oy Imaginary component of Yr 38b 
é Phase angle, defined in Appendix 5 A5—36 No dim. 
ys Phose angle, defined in Appendix 5 : AS-37 No dim. 
w Circular frequency of forced vibration ; A5—~6 1/sec 
i Inaginary part of, w=aP,H2/1%R2 37b L/sec 
o, “:atical frequency, coteapondiag to 8= tia, 37b 1/sec 


and § = 8. 


=| pens 


\ SYMBOL. 
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SUBSCRIPTS AND OPERATORS 
REEERS TO 


Atmospheric, or ambient 

Critical 

Equivalent 

Ordering subscript (also isentropic, Appendix 5) 
Amplitude of perturbation, Appendix 5 

"ecess and recess boundary 

Supply 


Equilibrium condition’ 


Differentiation with reapect tot = wt, Appendix 5 
¥ fics ; —_v* ‘i 
Differentiation with respect to t ce t, Appendix 6 


Differentiation with respect to t, Appendix 2 
Differentiation with respect to z = ae text 
Differentiation with respect to x and r (or € = r/R), text 
Laplace transform of variable 

Dimensionless quantity 


Mean value of quantity 


Pleference ceantio 
sererence quantity 
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FIRST 


APPEARANCE 


5 
44h 
A2-1 
16, (A5—20) i 
AS~36 
6 


A2-1 
19a, 32a .& 
13, (37b) 
19 
A5-—6, 46-7 
1 
A6-7 
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APPFNDIX 1 


THE PRESSURE INTEGRAL, 3.t 


In the integral 


R 
1 = { P(r) ede 
i, 


ip which 
rt 
ite) = pe (1 ~ Floren) 
and 
2 
E - P.. 
R, 
let 


2 
Pro 
2.2 a2 
R Peo Pat bares 
r=R cox e 2 
FUR 
: 2 
P 
fo 
peop? 2 
R\ fro ‘at WF 
do" ho e° 
t 
VR, 
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(Al-1) 


(Al~2) 


({A1-3) 


(Al-4) 


\The method employed here iu derivieg an expression for 2 is identical with that coutained in reference for? 
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so that: 
2p? Laid 
: : 2 2 
. Po-Pp % ARy é 
§=-P Pp? GF) * “(2 { Ge F a6 (Al~6) 
ro 
Rey 27 Jon, 
2 . 
Integration by parts and introduction of the factor pind then yields: 
wT 
z 
2P ; 
: aiaat = ; @(R) 
y fe “Tat? 2 62 _ 2 2,2]. 
a2 Pir? {| oye ae Ss Pe a6 (Al-7) 
Fore 32 R, v7 v0 
a(R,) 
in which 
; ay 
é(R,) -(3) 
Fc 
(A1-8) 


p % 
G(R) = -= (4) 
Flo Ff}. 


Values of the probability integral 


are listed in standard inathematical tables. The integral can be used in determicing the load 
capacity, I., of the bearing: 


2 
PR 
1. = nR*p ( 23 aaa se +1) (A1-9) 


A good approximation of the load capacity, L, may also be obtained in assuming a linear 
pressure profile. (See Figure Al-1.) 
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FIG. Al-/__ PRESSURE PROFILE AND COMPARISON WITH STRAIGHT LINE APPROXIMATION .— 
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APPENDIX 2 
SIMPLIFIED STARILITY ANALYSES 


In the simplified approach to the stability problem [2,5], the flow in the slot is treated in 
terms of lumped parameters. and ihe dynamic pressure profile at any instant of time is assumed to 
be similar to that which would result if the bearing load were applied statically. In the case of 
the circular bearing, a further simplification is made by assuming a linear profile along the annu- 


lus. (See Figure Al—1, Appendix 1.) 


‘The perturbation equations, which correspond to equations (13), (14) and (15) of the more 
exact analysis, can then be expzessed in terms of the variaticn of gap width, h. and of the recess 


pressure. p_. as follows. 


, oe o 
mh=Ajp.+ JJ pad +6) = Alp +£t0 (A2=1) 
A 
a 


Ah + App, =~ Agh ~ (Ag + Ag)p, (A2—2) 


F hy 
cin whieh 


A. = area over recess 
A. = area over annulus 


A, = “equivalent” area 


wee le tep hen, Vr 7 
AS =x 11 rata 1) ———___"__ (A2—3a) 
aX, 


A, : = mn? I(e) + (F) ¢ | es 
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while the coefficients »; are as follows: 
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: : (A2~4b) 
- pt, We Pe : : | 
al (sr “op ty 
ee “ro 


5 ap 
he Po = Par 
Combining equations (A2~1) and (A2~2) we obtain: 
A, tA AL ALY AL A 
aos 3 5\ + . 
re P(e ee (A2-5) 
AR m Ag m dy, 

The characteriatic equation jaa cuhic- | 
3 2 aock ; 4 
s* +A,3* +A,st+A, = 0 (A2-6) i 
' 
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and scakility conditions require that all coefficients A; he positive, and that 
AyAy > Aj. “(A2-7) 
Vhe inequality (A2-7) reduces to: 
Ay. A 
1 4 
> z (A2—8) 
Has AG ee 


from which the critical recess depth, indicated with broken lines in Figures 3 and 5, was calcu- 
lated. : 


It is to be uoted that the coefficients A; are analogous to the coefficients C, and K; in 
the more exact analysis, but that there is no one-to-one correspondence, except for A, and Aj. 


Next, consider the relative, axial motion of two thrust plates, one integral with a rotor, of 
of mass m, the other with a support, of mass M, resting on isolators, or on an elas.'= foundation. 
We associate with the supporting structure the lumped elastic and damping constants, K and c. 
For the resulting, two-degrees-of-freedom system, shown diagrammatically in the appended Figure 
A2—1, the following set of equations may be written: 


mx, ~A.p, = f(v) (A2—9) 
Wx, tex, +Kx, + Alp, = 0 * (A2=10) 
ALB +Ayp, tAgb+CAs +Ag)p, = 0 re (A211) 
in which 
h = Xy-x, (A212) 


Applying the Laplace transformation to equations (A2—9) through (A2—11), and assuming initial 
rest conditions, we obtain: ‘ 


A = 
oe yy Ree ce f(s) , 
x, m= F--: (AZ=13) 
m m 
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oe : \, : aha ; 
(s - ss — a) yy _— po (A2-14) 
\ : , 


in which 


The characteristic equation is a quintic: 
Ags® + Ayst + Ags? + Aya? + Aps tA, =O (A217) 


in which the coefficients, A;, are as follows: 


Ng = Ay 


CA, + ( Ag t Ag) 


- 
"t 


(ord Mo + COA, + As) + (B + b)A, (A2—18) 


> 
wo 
" 


> 
nN 
" 


M2 (A, + Ag) + CbdA, + (B+ BIA, 


A. = h(a, +02.) 
ry a a 


oer 


= 2 
A, = bOPA, 
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Nithout loss of generality, one may assume that 4, > 0. The necessary and sufficient 


conditions that all roots of the quintic have negative real parts and the system be stable can be 
stated as follows as? : 


! Ay ~ 4 (A2~19) 
: \, As! ai 
a Ns As (AA, > Ag \,) > 0 ‘ (A2~20) 
ji, As e | 
ee NA Ns ~ ASCASA, ~ AQ AL) + Ag, > 0 (42--21) 
\\ \g} 
\y A. 0 0 
oe A, 4, A 35 ” 9 
< aes =~ [(AgA, ~ Ag Ay)? + (AgA, ~ Ay Ay dy ] > 0 
ea (A2—22) 
Se = AS Ay > 0 {A2-23) 
Since As. A, and A, are generally positive, it is sufficient to examine the signs of 4,, 
oy and oy : 


Two points in the latter part of this analysis arc of particular interest: 


a. If the damping coefficient. c, is zero, the stability conditions of this system are identi- 
cal with those given in the inequality (A2-8) for the case of a rigidly a 
bearing. 


b. Since A, = A, is the only cvefficient contet sing ihe recess depth, 5, as a parameter, 
the determinants A, and A, yield a quadratic and a cubic in 5, Consequently, unlike 
in the case of the Gndanped: rigidly supported bearing, more than one value of 3. may 
correspond to the same value of the pressure parameter, P,~P, bPig 


In the toregoing analysis, it was assumed that the external dampiny for tke rotor, ot mass 
m, was negligible. This, indeed, would be the case if gas bearings were employed as a meuns 
of lateral suppor:, but co Uiffivuli;, vines chaz at ulpetcaic one, is involved in adding a damping 
term, —cyk,, tv the left hand side of equation (A2-9). 
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APPENDIX 2 


pert avery an crawernrne 
eee ee tte Nae APG Sore ate 


Soar gs — a a ac aE | ‘ 


Des ce 14.7 psia 
P. = 88.2 psi 


m7, 4.65 x 1073 [b-sec? Gu.‘unit length. 
‘1 = 530° F aba ; 
F265 «1079 tbhesec ‘in? 


%* Feve = hd 
v 


2X, = O.2in; ay = 7.54 1074 in? 
NW. 2X = Sins ay = 7.54 = 10°74 in? 
I. 24, 7 OL in; @y = 6.86 = 1079 in? 


' m5 | 2.21 
' t3.0 | 2.95 
i RO | 5.64 
{ s3.2 | 9.70 
i 85.0 | 15.0 


~y 
i) 
uu 


fo] 
> 
Qo 


ae 
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APPFNDIN 4. 


CHRCTE A GhUME LAY - DETRRMENATION OF CREFICAL 
RECESS DEPTHS. - TABULATION OF CALCULATED AND 
FXNPERIMENTAL RESULTS. 


Nhen explicit expression of the coefficients (;; in equation (16b), are substituted into the 
chacacteristic determinant (315) and the complex variable, s, is taken as equal ta ia, the latter 
Cad be written as: 


Hy =O (As-1) : 
in which 
ee ee 
hy joe - — 2 et) 
Pn? ttt? RS 
toor t 
hry O Fy Fie 
Ray 
(A4~2) 
ee \ alo 2(F)(2) 2a ; 
ay 7 2 GD $2 [J] —] fio 
} 4 \p R} aé | | | Ss 
t i 
2 2 
Pat R * a) Pat / 
-2})+{—~- _ = (LD +]b-[— ! 
Rol WR] o€ Po fs | 


epee 


| ci 
Gay - PO Rimes {o.bp = 2p — — U1 “ia B)- j A4—2 
rege Ga) ae wi 
\ ath i , . ; 
Me Ge (A Zt . ; 
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{an fPaV rw! an 


and 


TL ey Lae 
ay [tajery + eytiag + Heyey O04) 
Sete, ee Ne ee a ee 
2 2 
Oy + My 
€=5 


HO a, (Sy Sey + Uy Oy + by Oy — SM) 
as¥ = tS HS 


ot 2.2 
Gy + Us 
I if f=5 
wn Oy id 
a 
is 
7 oy +e 
Wl Wl Fas 
In the foregoing expressions: 
12pR ow 
a= 
Pou? ) 
ro 9 
2 
-(#) | 
R= Pro (A4—4) 
lage oe 
of, Tf 
t 


vo 
7] 


R 
f Pi(r)rdr (Appendix 1.) 
R 


r 


and Gig}, Wigs are the variables in equation (405), solutions of which, for discrete values of a 
atid for (3 = 0.596481, are tabulated in Appendix 7. 
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Cet the bearing geometry and the rarameters Pay Poe Tor ¥ aad p be specified. The gap 
width. HU. can then he exnreaned in terms of P.. the supoly preanure. by equating the steady- 


state muss-flow rates throngh the nozzle and through the slot, e.g.: 


9 | as 
ripe -1 P. | 
' a ad 3 2 Y a 
Kee mas un? ~ {22s aa) ni(e-) (A4—8) 
i2,: T, 4 ‘T, yl Ps 


Given a, the unknown parameters which appeur explicitly in the characteristi: determinant 
are PY and 6, which must be selected in suck a manner as to reduce simultaneoualy to zero the 


real and imaginary parts of the characteristic equation. (e.g., §, + & = 5, in equation $3b,) 


Tuis reduction was accomplished by trial and error, ae shown in the following sample 
calculation. 


Sample Calculation 


Given: 

m= 6.477 » 107? Ib-sec /in 

Bs ~ O5in 

R = 2.5 ia 

Pi, = iT psia 

P., © 73.5 psia 

TS + 530° Fabs. 

£2 2.47 x 10° in?/sec2 PF 

yoo 4 

> 2.65% 107 Ibesec/in? | 
Tet 

a = 0.950 | 
and assume | 

a ‘ ‘ 
Y — 79.6 peia 


THE-FR ANKI iN INSTITUTE © Laboratories for Research and Development 


Then 


hes adel, 
i= O12 108 
and from Appendix 7 have: 


ty an OLSTaS 


OS) = HARV INN. ZR) - 06RET 


tay 


From equation (A4—3) have: 


aH 
2— (21) - 0.26404 4 0.2090 
ag 
ee) 
2— (1.5) = 0.4860 + i 0.3788 
dé 
Hn) 


2 — (1.1) = 0.5312 -i 0.1378 
a& ; 


Faquation (A4—5) yields: 


P 2 
\ = oat «0-4 : -1) 
° P 

ro 


(5) = 0.915) : eG) i 


0.8906 


= 0.6019 x 1075 


fo taat 
H” = 0.4096 x 107? w, = 0.2465 x 107® 
. He = 0.1825 x 10-5 
H, = 0.1351 x 107? 
and 
fall.) = 0.3082 x 10712 
78 
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i Dette ee OL RAD < Shae 


I-A2049-12 
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The heuie of Jie determinant (\4—1) then become: 
\ ty = § OS5TT 
| Gry > 9.860 > 167! ~ § 0.5983 
| Ty, = 0.8620 x 107! - § 0.4060 x 107! 
Bay > - 0.3295 x 10°) - 5.03571 x 107! ~ 5 0.2644 x 10? 8 
Gy, + Oleots 10h - 4 G.2020 
] Try > 0.4550 > 10! | 6.9377 ~1 0.2644 x 1075 
Fy © 0.5000 » 107! 
(4 = 0.7600 » 107! 240.5101 
| Try ~ 0.7600» 107! 4 § 0.5601 
Fy, ~ 10.5000 y 107! 
Gy, = 9.1812 » 10h 4 § 0.2143 
Yo, = 0.1812 « 16" ~ 1 0.2643 
For a = 0.050 and P. ~ 79.6, the characteristic determinant (44~1) reduces to: : 
WM = TT = Thy Tey 
- alS-¢,) + ib(5 — g,) } 
t 
= = 1.073(5 - 0.3668) + * 2.979(8 — 0.3639) 
} 
t 
76 
{ 
t 
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Hiesults of similar calculations for other values of Pare as follows: 


« . &) 
BLO ORK" 
80.0 0.3274 
79.8 0.3455 
79.6 \ \ 0.3668 


A plot of 2) amt ga vs. Pl, for @ = 0.050, is shown in the appended Figure A4--1. 


Be 


0.3724 
0.3675 
6.3655 
11.3420 
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0 380 
0.370 
a mf: ee es a4 
cert a dE 
pop Pep ek fe] 
0.350 
0.340 
2 oe: 
= 0.5000 x 10 
ial 2 
; | Ba Ih/ Pol? © 
0.330 |-- B OG RURe 0.5965 
2 
m= 0.6477 210 EB= SEC 
IN. 
0.320 ——— 
79.6 79.8 — 80.0 
P, -SUPPLY PRESSURE (PSiA} 
FIG. AF-1. DETERMINATION OF THE CRITICAL RECESS DEPTH FOR a = 0.050 | 
| 
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1. n 

boo 

RAE Ayes 

en ee) 

oe) er ee te 

183.2) 0.332 

183.3 i 0.123 

i a4 0.135 

( 8t.2 t 0.146 | 

| i. om = 6.477 

ome pia 

'zo.4] 9.0803 

| 79.65 9.0930 

‘80.3, 0.0025 

ROK | 0.0966 © 

{ Si.0 - Gor 

| 82.2 9.118 

feogl oat 

| 83.1: 0.13) 

[826° 0.137 

1812: 0.146. 
85.2! 0.159 
87.2, 0.185 
89.6 0.219 . 
S18; 0.249 

| 93.81 0.276 


82.4 | 
88.45 


o1.9[ 9.228 


PRESENT 


0626 


[1] 0 
181.4) 0.108 


0.121 
0.203 


» 107? Ib-sec? ‘in 


Pe eras 


TABULATION 


ANALYSIS 


90,1 | Lot 
54.6 4 1.67 
$5.4 1 0 1.75 
1.6 | 1.79 
36.6 | 1.84 
31d | 1.69 


39.9 
13.8 1.62 
_9.60| 1.68 


{4.55 
110 | 1.72 
95.9 + 1.81 
a 0 | 1.85 
83.5 | 1.90. 
vay of 1.98 
49.0 2,24 
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* CALCULATED RESULTS 


= 6.055 in (ay = 0.237% 107" in”) 
ne - 14.7 psia 


RP - 728 nsig 


T= 830° Fabs. 
= 265x107? Ibesec ‘in™ 


CT ne ow ee ee a 


St 
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TARUL ATION OF FNPERINENTAL RESULTS 


Twelve 0.0167 in diam. nozzles (Total throat area = 0.263% 1072 in) 


T+ W21b (theoretical load at P= 5 PL, = 73.5 psia) 


loom = €.477 » 107° th-sec?. in; = TL74+ psia; T 537° Fabs. 


" 


’ 
Pate o 


536” kabs. 


o 


Nae a aaa ee aa aera: | 


11.98 » 107? Ibesec? in: PL, = 14.63 psia; 1 


~ 
aaa 


| 
{ 
| 
i 
} 


itz) 


| 2297 
2214 
2185 

| 2153 

| 9126 

1 2076 

{| 1987 

i «186 

1802 

227 

1689 

1630 

1487 

1424 

1380 

1536 

1307 

1277 

1232 

1334 


_ | 1346 


tf) 630 i Ay 
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APPENDIX 5 


YEE OPP RGGi CHARCESG AND DISCHARGING 
OF A VARTARBLE-VOLEME RECESS 


The purpose of this appendix is tc invesstizate, at least quatitatively, the adequacy vi the 
asaunption that the process of alternate charging and discharging cf the variable-volume recess 
‘unid be represented by the isentropic pressure-density relation P pY = constant to a reasonable 
degree of anprosimation. 


Tn the course of experimentation, small amplitude auto-oscijiations and sine waves have 
been viewed sinultaneousts and compared on the display of the cathode ray oscilloscope. Ob- 
served differences in wave form were slight. We shall consider, therefore, the case of prescribed, 
sinusoidal motion of one bearing surface and. making a series of simplifying assumptions. similar 

7 


to those in Appendix 2. utilize the energy equation, instead of resorting to an a ba assumption 
uf the pressure-density relation in the recess region. 


In what follows. pertinent equations will be deve loped and linearized. Resulting perturba- 
tion ainpliviles of pressure and density, as well as phase angles, will be related to those obtained 
when Pp? = constant. Numerical examples will be given for parameter values corresponding to 
representative points on ine double valued, experimental stability locus I, Figure 55. 


Tt will be assumed that perfect mixing of the entering gas occurs. but that fluid velocities 
in the region V_ are negligible and properties uniform throughout at any instent uf time. Adiaba- 
ticity and eanali of temperatures T_ and T, of initially steady flow will also be assumed (e.g.. 
Tle, = 1). The flow rates and their perturbations will be expressed identically as in Appen- 
dix a ‘ioplsing the quasi-steady flow and lumped parameter approach to the problem. The fluid 
will be treated as a perfect gas with constant specific heats, e.g.: 


P sees ¥ 
— = *T; #C -c: ot = 
; J ¢. = y 
so that 
1 
ue ce T = Z 
= y-1l ? 


rt 


where is the internal energy and h the enthalpy. 


al 


fs 
i 
I 
i 
t 


~ 
S 
Pon 
iy 
| 
is 
Sg 
t 
| 
‘ 
' 
t 
| 
1 
i 
| 
dq 
a 
J 
\ 
\ i 
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Ibe continuity and energy equations for the system shown in Figure A5--] are as follows: 


Wie ee es eS (AS=1) 
run de de Pr 
Woh toh e 4 vu) P. (A5—2) 
l=s 2-2 at Pride + he 
which a be written as: 
dv, , ape 
ae YS (A5=3) 
P dV, iP, 
Rota ieee ee (A5—4) 


a oar as 2a 


The following dimensionless perturbation quantities are now introduced: 


, = Wo + Ww) = Wa +w;) + 
wf, = WN + Wo = Wi 4%, 
a a pf +p) 
Pe = Pr + 8 = p_ (14a) 
Wo= Mo +h = Hol +h) 


r 


t woke 
Lf) 


Let also h be prescribed as: 


h = h sinet = H esinat: = - <1 
o a 1 
o 
Su tnat 
& - € sunt 
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V = AXS+I) + Ake Vo + Ah © Vi + Ab); B= 


_ (A5=5) 
Il, 
j 


f (A5—6), 
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‘ , Also. as in Appendix 2.. we take; 
< ae 
: pea beol 
a¥ I ) Wy . : 1 Pro 
| oP, Pro 8 OL By = 
fo 
i, ber (<*2) WL e 2 
As -~J ~'=—-h, dg = ee 
5 \p, ip Po 5 - (Es) 
Pro (45-7) 
ee a ws \ 
={- = od Ay= 3 
’ e Fae | ee ae a \ 


\ \ 
A, t+ Ag = — (X, +X5) = —A 
P ; 
ro ro 
; so that with: 
W,-%, = ~Ap-4,h (AS-8) 
equations (45-3) and (A5—4) yield: 
\ a+ Qe - , P-Op = 0 . (AS—9) 
. + QAp = - Qa, esint - Pecost (A5—10) 
in which 
v 
: a ae 
Pro ro ~ me | 
{ (45-11) \ 
: i 
B 2 No } 
6 + Hy H 


d 


qo ae ah denotes differentiation with respect to dimensionless time, t. 
t ao ait 
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| | | 1-A2049-12. 
‘S \ tn operational form, equations {45-9} and (AS—10) are: 
\ : . as-(Lp o)é £6) . (5-12) . 
: Ae (~p+o Pp _ (A5- : 
Da + QNp = -e(AN+OK,)sint (A513) Og 
where 
d 
n 5 
On the other Land, an assumed relation 
wt a P y : ~ 
ao Peo fee (AS~14) 
; ro Pro ‘ 
vields 
tay | (A5=15) 
a 
so that equation (A5—3), with the aid of equations (45-8) and (A5—15), reduces to: 
. + QAP = -Qi,esint - Be cost (AS-16) on - c 
ac ’ ‘ 
oF: te | 
( : b+Qh)s = -e(pp + Qh,) sine (AS5~17) ; « 
a ‘ | ; 
Operatiot: on equations (A5—12) and (AS—13) with (D) and (D + Q), followed by subtraction, and Og 
operation with (GA) ard (1.’y D + Q), followed by addition, yields the following set of equations: pe: 
pL Ar? +008 +5) + 07%, Je ot | 
B= 9 | 
[ P0144) 0+ | : | 
; 
} ~ [(Q?X, -A) + i0(8 +%,)] tg 
4 Sy eos { 
te = {————____—__——- (cost + isint ) "(AS—18) 1 


(o?i- 1) +i9a +A) 
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the eer band. for cometian A037) we antaine 
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(A519) 


i (i emitted A Ny ict (cost + isint) (135-20) 
(20-07) eae Bee 
y : y 

In the forecoing equations the subscript i denotes the isentropic relation assumed in ohtaining 


eration LG=1T, and. cefers to the imaginary part of a complex quantity. 


The steady-state pressure and density responses to a sinusoidal displacement input: 


ho: ¢sint 


ate as follows: 


a ae he ae we (aes 3 
: bE) Ay PG “AD AU way ar y- )| 7 
Pr A cs (LISS 


f = ( a 7 Dy ~ a 
«[lQra, - #2) OFA - ra + (B+ AMI + a) | 


Ce sit 
wD. 


iy 4) +d. a? 


{45-215 


aac ee — < oe Sine 


meee 
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atl x 


note, once again, the form of the parameters { and (): 


H 
R= = 
+H, 
and 
WV 
os — 7 
Pro Vio 


T-42049-12 


(A5—-22) 


(45-23) 


(A5=21) 


The expressions for Pf and o contain the parameters ¢ y, Ng A, B and (), of which ¢, + 
, ate constants and \ is fixed for any given set of pressure ratios PP and PP. We 


Prior to solving equations (45—2]) through (A5—24) for a set of numerical valu2s, it is of 


interest to consider the asymptotic behavior of f and & for limiting values of the ‘che rge” 


parameter Q, We shall consider the cases when: 
(a) Q) is a very large nun ber 
and = (b)_ is a very sinall number 


whic Bo romaine finite. 


RA 


Case 


‘ 


so that 


but since the assumption (PP) 


aw es -¢ 
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th Eis wees Lone 
I} suations (VG=21) and (A5—22) reduce‘to: 
3 Ay le j# Sie 
! por ne re sint (AG—25} 
( 3 
: 7 
woe me ee sint (A5=26) 
A 
4 
as P a - 
Pe ee i eg te (45-27) 
a o Pro Y 
We note that equation (45-23) also reduces to: 
ay ad 
P= re 5% sint (A5—28) 


(p,’p;,)” prescribes the p, and G, relation, we have: 


Pp, P. .P a i 
2 =y or tan fe vero = ay (A5~20) 
9; % Pro 


In case (a). the following conclusions can be drawn: 


I. The pressure-density relation is nearly isothermal. 


2. The pressure and density are in phase and both we 180° out of phase with the dis- 


placement h = eninT. (e.g., increasing with decreasing gap'width.) The coefficient 


is a measure of the “‘statsc stiffness’. The phyzical significance of the foregoing in 
that, if the motion is very slow (quasi-static), or if the mass flow rate is relatively 


large and the mace rortent of the recess region is relatively small, the changes vtcus 
"Nearly isothermal ly. , 


a7 


a ee ee 
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tose thir thas veew seal. 


Porcatiens CUS=07Y and (AA=22) reduce to: 


p= ~eyjisint : (A5—30) 
@ = -efsint (45-31) 
so that. 
p ~p "3 
Ve aiah; Ngee whe: ae (45~39) 
‘ Pea be : oa 
I quation ((5=23) alse reduces to: 
PB, = -eyfRsint (45-33) 


, p 7P. 
2 ey. ee se ee (A5~34) 
% 9 Pro 


which vas the assumed reiation. 
In case (bh), the following conclusions can be drawn: 
1. The pressure-density relation is nearly isentropic. 


“2. The pressure and density are in phase and both are 180° out of phase with the dis- 


placement h = esin f. (e.g., increasing with decreasing gap width:) The coefficient 


> 


is indicative of the “dynamic stiffness’’, such as would be the case of gas in a closed 
cylinder and piston system. undergoing isentropic compression or expansion. The 
physical significance of the Zoregoing is that, if the motion is very fast, or if the mass 
flow rate is relatively small and the recess mass content is relatively large, the changes 
occur nearly iscntropically. 


ee een oH se 
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To correlate the analysis with experimental reaulta, the numerical examples were selected 
iv correspond aaproximate!y to two points of the experimental stability locus I, Figure 55, for 
which tne data are as follows: 


Sa me a eA fois 
2 


pews oey in ee eee ma a se 


i PE Pia | 915 
Dre | psia 75.9 
‘ ! i 147 H 
Pat \ psia 24.7 | 
an babs, | 537 | 
Pro | Ib-sec” ‘in’ 0.5722 x 1076 
WV, ! lb-sec ‘in 10 x 1076 
HN, | in | 1.64 x 1073 
A, in@ 0.7854 
= 35) 
i, | : 1.933 (A8~35) 
Wea a | 2.078 
ae i { 
t = | “4.011 
ry - | 3.000 
- oy 2 1.4 (1/y = 0.7143) 
1072 
{ @ | rad/‘sec 
AE 11.50 x 1073 217.3 x 1073 
| | | 0.7153 6.204 x 1072 
| 6 | | 0.1120 5.928 x 1073 
ee Se ay ee: el CE NY | 
| 
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; ae 
Pyare ft 


rn 


oh sin (tT~é) ; 


tania + 


Om 


Pm 


cf 0.792 sint + O62 tcost) 


0.1624 


— - 0.2361 
~ 0.6792 


186°. 13.45° = 166.6° 


0.985 


ee sin(t 0) = ef 0.5291 sind + 0.232cens7) 


m 


tan OF 


ee |: 
~ 0,529) 


Ine ~ 22.87% = 157.1% 


> 050K 


eeu = OF 


> efi, sin(t -¢,) - e(-0.7135sint - 0.1386 cost) 


0.1386 


opr. = 0.1943 


180°- 11° = 169° 


0.7269 


! 
| 
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(AR) 


(45~37) 


(\5—38) 


(45~39) 


(A5~40) 


(A5—41) 


aa 


4 \ at 7 ‘ _, : ro Neg : 
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Teampe 2. 
Pp - ef. singt 8) ~ e(-O.091L sint + 0.2218 cosT) 
; 0,22 
tai go = eral e. ~ 2.468 
| aa ata (A5—42) 
| | G ~ 100° 67.93" =. 121° | 
| ; FR, » G2H5 | 
| 5 = 8 snit-u) = e(-0.0611 sin t + 0.1625 cos t) i 
eG . ee ee 
| ; - 0.6611 (A5—13) 
a ee - “= 100°-69.4° = 110.6° 
a 5, = 0.1765 
; o-w = 1.5° (A5—-44) ‘ 
F, = eff, sin(T-4,) = (~0,0923 sin t + 0.2322 cost) 
tig, MO DARES so laps 
land; = 9003 2.516 
é, = 180°- 68.32 = UL? (45—45) 
f,, = 0.2399 
iw Sp, 6 OT85 
wl y Pri vie J 
¢,-¢ = -0.4° (A5—46) | 
d,-u = 11° (A5~47) 
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; ( wuparison of resaits of examples | end 2 shows that pressure-density phasé engles are 
+ + ‘ “ -o . ‘ + . 
— relatively small AP pean Hey tole = nk Too respectively, When the assumed relation is 
. 4 3 . ; < 
i ai Pet a J tha: ie eenivatent to prececihine that the densite and pressurc he in phase, 


’ F rtd eae 3 
Ti cane ple Lo the ancle between pressures is approxisnately 2.4” and that between densities 11.9". 
é o z, 
he exanple 2. the corresponding angles are very small; 0. f° and 1.1" respectively. 


Since all uncle: are relatively small. so that all densities and pressures ure “‘uimost in 
i phase’. a comparison af anplitudes is of interest as a criterion of validity of the isentropic 


rreren 


Promevareple Powe have: 


Pri ‘ Oni . 
=. 1.0% and — = 0.904 . (A5—48) ~ 
P, on 
alse: 
Bi . Pri e 
——. = 1.2. as compared with —— = y = 1.40 ((5—19) 
at ni 
re 
In example 2 the ratios are: 
Pai Oni i 
=—— = 1.02 and = = 1.01 (A5=50) ae 
Pm om , ; 
alse: ; : ‘ _ 
Pa Pri . iss 
— = 1.37, as compared with —— = y = 1.40 (A5~5)) 
Om Omi 
In the foregoing examples, all parametere were identical, with the execption of the fre- 
wueacies and the recess volumes. The ratio of the former was of the order 2; that of the latter i 
of the order 1 20. Significant!y. as might have been anticipated from physical considerations, } 
the pressure and density amplitudes of example 2 were burely 13 of those in example 1. 
: \ 
ats ee : r aes 
For a sufficiently large volume, pressure and density changes in the recess, associated : 
f with small but finite motion, could become negligible. A very different condition, of course. would © a 
: prevail in the bearing gap. due to the combined squeezing and trapping of the fluid between the | 
nerrowly spaced houndaries, . i 
fs ! 
‘ H 
fies tits of calculations of examples 1 and 2 are illustrated in Figure 45-2. ; 
i 
4 Ne 
iy ‘, 
m i 
i} 
\ 
ea = An a %, 
‘ ; i : : rat cart ncn meetin FY 
’ ib ay Sed 4 = \ Re, 
hee 1 . Py me S 5a oe iets 
? x 


=: wo 
J oe a ene 
ve ia eer 7 a wok Pome - - ae ‘ 
com Se eee ta ee T+ ++ 5 » -| f . i 
- 1 
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ae FIG._A5S-1 SCHEMATIC DIAGRAM OF VARIABLE- 
VOLUME BEARING RECESS 
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SPI ee a 
— = 166.6" a 
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7M 2157.18 
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4 £ ; ‘ ¥e a, 7 €im : . q 5 a a : . 
<3 i ; P 4 . : oe : ; 
. € Elin 
See EXAMPLE |. | 
~EOmi 
: —€P mi } ; 
i | | 
EXAMPLE 2. 


“FIG. A5-2. _AMPLITUDES AND PHASE RELATIONS OF OF PRESSURE AND DENSITY . 
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The purpose of this appendix is to show that the differential equation (4) of the main text 
pres des an adequate description of fluid-dynamical phenomena in the bearing slot. provided orders 


ct oaeenitide of characteristic parameters fulfill certain requirements. The development follows 


sosely the method of Flred [15] vho employed a small-parameter technique in the derivation of 


catitions for gas lubrreation. 


Ht will be assumed that the fluid obeys the perfect gas law and that the specific heats, 
Viscosity and thermal conductivity are constant. Continuity and constancy of temperature of the 
boundaries will also be assumed. In the range of presaures and temperatures considered, the fore- 
#0ing assumptions ace reasonable for gas flow in a sarrow slot between metallic members. 


The momentum, continuity and energy equations, stated in cylindrical coordinates, for 


avail; sxmmetzic. laminar motion, and the equation of state are [65 1: 


ree 


at or dz 


dr, 


- lid 


3 


a (au 
ar \dr 


r 


+ 


7, dw 
cae 0 


P (2 ia 
dr? r 


or. 


w= 


(46-1) 


) (46-2) 


(A6—3) 


(464) 


(A6—5) 
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c 
‘ 5 : ) 
in which alk > oe) ee oy - 2 are constant. 
Cy : 


Livi teenie Aes hiner ve 

Woz 4 u w: 0 
yoy 

AL 2 Wr); 1 Q 


(A6-6) 


Denosing reference quantities with( )* and dimensionless quantities with (7), the follow- 


ing set of variabies is forn ed: 


Moon ptr* = canatint 
o> pp 
I vot 
9 
hte 
p- Gene 
PAN" 


Wo. Uf), whee Midis oO. 


ro- LT, where [, «= R-R, 


(16—7) 
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To the foregning «+ 


the dynamic lulvication equations, 


Toi 


ve 
u 


rs is thc small parameter used in developing 


Ms, 
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In transforming the equations of fluid :wotion inte their dimensionless equivalents, the 
following eelations between derivatives hold: 


at ? 

a ] : 

oF LL F 

or I. 

a7 fn* 

50 hat: 

9 ~~ DO 
or 1. at 
eee Oe 


in which (°) = * 2 2 . 


= 0 


ve f 
aoe Ae : 
4 


aro? 
Cae RO re 
* a 2 


a 
dz 


“ab 


Qs ja Qiw 
sig =| 


Uy 


eee eres ee eee 


(A6—8) 


F 
eee 3 i a ait ge a = 
f ae 2 . oo ‘ paseo rT ae mn ee tad a “ as 
- enue + 
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The dimensionless ‘oun of equations (46-1) through (A6—5) then is: 
oo = y = ey - +t I toa P 10 fat nv lL aw 
aa - wo Z - uu cress Tar? — — tert —- = 
w NE : Ms oF eo dF 3 (ae F OF OF 
(eo i ae TF lt L we ng 
Vee roar 7 Pe {< ay; } 
: 4 = 
18 Q PLN AY -aw il : ] av 
Vogt Jae PL Az ei | A f AF 
. eta ‘ams _— 39 om ~ 6. Fs 
bbe (8.1 aR). (EE LLL LEE) er0 
3 f oFsor FF * a or- Ot OOF? dz" 
1 
me fe _fs oe ,1 orgy. 1 © (pH) - 0 (4611) 
rs 2) a f 3% 
a Rata {-o\er av 1_1| oP. fleet s\eP aP 
ee At Le fz) st .o A -L a, (te-fz)e ict 
ite lot f f OF oF | e-l at f f dz OF 
: 2° - ; ~~") 
ap a RL Os Ge Pe a 
BC RSL OES OTe ek 42 ag? 


2 ey . 9 ; 
1 (3) f? \ dz 
‘Lhe boundary conditions are: 
At z- Q : fi - W 


At Z- | ’ 


Uae al 
P 
- 0 


(A6—13) 


(A6—14) 
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Consider now solutions of equations (AS—} through (AG PQ) which can be expanded in 


‘ ins 
Tevmes ofthe small narameter¢ ao follows: 


{ ,o ee 
j ae co & 
vo that: 
tow 
os i ‘ 
Mae” ae, ap 
Po oT 


ok) aaa 


Wy + . 
HedW eee 
eV Ey wees 
Sc incn cvs is 


5 nee geld Se 
P*| 6, 6, 


1 


| 
| | (A615) . 
mei 

an 


The functions appearing in equation (A6—15) depend in general on t, t and 7 and are in- 
dependent of... ‘The zeroth order functions are made to assume the boundary values (A6—14), so 
that approximating feactions of higher order vanish at the boundaries. e.g.: 


sy 
4 
Q 


At 


NI 
1 


At 


are: . 


l. 
v 

"0 = 
6, = 
» 

Wo = 
9% = 


~e 


" 


O07, 1 
~— 
or f 


J 


(A6—16) 


Substituting the series (AS—15) in equations (A6—9) through (A6—12), the zeroth equations 


(A6—17) 


vr. . tees + - . . t-5, + . - 
‘tower ute dittering by two are taken since exponents ot ¢ in tne differential eysations differ sv twa, 
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(A6—18) 


= 7, (7. t) only, the continuity equation (A6~11) 


wpinties CAG6-181 that 7 


(A6—19) 


(A6—20) 


where in the penultimate step, use was made of equation (A6—17). 


Integration of equation (46-20) with respect to Z gives: 


72 2Q=2 OF i ey 
Me gy HE SEs oeele-c) G(R T) 
2 2 ot 


und. applying the boundary conditions (46-16): 


2 OF, 
Le. s ae gripes La (A6—29 
2 at 


av that equation (A6—2:) becomes: 


ec *2 = 
a A{ Se a eee Lee (A6~23) 
2 Bs at 
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be the reference velocity. corresponding to the steady-fliew condition when f = 1, equation (A6—23) 
can be writtes in the following furai : ‘ ~ 


| Ps ah <4 Rie Te \2{ nf 
fo4 Piva ae bea _Z) (=. me nf (A698) 
ge OA Ee r/ at reyes 
in whieh 
Coy - * 
x; - Neate (A6~25) 
Vv y*T 
Sincss, genccall >. aa Np, !? ce 1. one may assume that to the first degree of. 
approximation , 
Tee = 1 (46-26) 


\ 


which is equivalent te the assumptios made in equation (2) of the main text, e.g., P/p =| constant. 
ft quatiog ‘\6—19) then reduces to: 


9% | 8 pee my. 1 aw, ; 
eee Ce le eee = 0 A6~27) 
or a a ae ar d 


tt {232 07% 
2 or 


+ ZA(F, t) + BF, t) (A6—-28) 


and application of boundary conditions (A6—16) gives: 


. of? a , ; 
Nee 48 het | (A6~29) 


101 


ene 5 


my: 


THE FRANKLIN INSTITUTE e Laboratories for Research and Development 


F-A2049-12 
wo that: 
| TS = 2 et-ae? (A6~00) 
on eet ead 
Subatitutinn of eqnation (A630) in equation (AG—27) vives 
nF , eae n= NV cé 
Se gees fe ae ee! fe ee a (A6--31) 
co A 5 | ° OF f ° az 
and integrativa between the limits Z- 0 and Z =1 yields: 
po eee 
j OF, 
£2 (re, St) - 12 2 (Fm,) (A632) 
t oof oF 


aa | 


which is the dimensionless equivalent of equation (4) of the main text, whea it ia aseumed that 


* is a goad approximation of P 

The error involved in neglecting higher order terms of the expansions cauld be estimated 
only if zeroth solutions were available; through substitution in the initial equations, for example. 
It ia self-evident that the degree to which zeroth terms of ¢-expansions approximate the actual 
functiona depende on magnitudes of other coefficients which appear in the dimensionless equa- 
tions of fluid motion. ( part.cular importance is the magnitude of f which characterizes the 
effect of boundary motion. 


It was tacitly assumed in taking Taylor’s equation {1 jam the starting point of the analysis, 
and subsequently implied in the linearization thereof, that f was arbitrarily small, In doing s0_ 
one was encouraged by the fact that equation (4) of the main text is exact in the limiting case, 
whenever: 


. 2 
L— «, with f = a Sf tiaite 
v dt 


or 


* * 2 
H*-+0, with a (4) finite. 
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APPENDIX 7 
TAREE VEION OF COMPETER SOLUTIONS 


‘This appendis contains a tabulation of con.puter solutions of the set of simultaneous 
equations (ii) suige:t to the boundary condttion C415), ‘the solutions are for a fixed value of 
A OSOGNRE and for discrete values of the parameter a. The ranze of the independent variable 
is ti €2 Sand the interval size is A Z- 01. The numesical method employed was that uf 

Hw ae q © OBasince and Dimsdale and is contained in ‘Numerical Analysis”, Kaiser 8, Kunz, McGraw-Hill 

Rook Co., F257 np 26-208, 

‘The Cillowine code is to be used in reading the tables, ‘The first digit is zero if the 

‘ number is positive and is a minus sign if the number is negative. The second and third digits 

stand for 50. n, where o is the exponent in 10". ‘The remaining digits represent a decimal, 
which must be multiplied by 10" in order to read the required number. 
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FIG_A7-1. PLOT OF COMPUTER SOLUTIONS OF EQUATIONS 406 FOR a = 0.093 
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